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ABSTRACT 


The main aim of the study was to develop a tephnique for 
obtaining large-scale, synoptic views of ecological variables 
such as plankton blooms which could be used as a convenient 
operational tool for modelling the Baltic ecosystem. 

Using all data from both the Stockholm archipelago and the off- 
shore test site in the Northern Baltic proper and, from LAKDSAT-2 
tapes and ground truth measurements, a strong correlation was 
found between the product of bands 5 and 6 (5x6) and chlorophyll 

The bluegreen bloom of Nodularia spumigena , typical for Baltic 
conditions was traced on two overpasses both by liANDSAT-2 and 
by use of helicopter tracks. 

Indications of current gyres, upwelling sites and discharge 
plumes were obtained from the "topography'.' of the blooms. Using 
the Swedish automatic instrmnent IRIS {Image Reading Instrument 
System) , Computer Compatible Tape (CCT) and programmes developed 
at the Sv/edish National Defence Research institute, quantitative 
estimations of chlorophyll a and Secchi disc readings were made 
and transferred to multicoloured maps using an ink-jet plotter. 

The combinations of techniques used should constitute a cheap, 
time-saving method of rough, large-scale quantifications of 
Baltic events. A continued processing of collected data will 
result in- a rude quantification of a total bluegreen ' bloom 
in the Baltic, Its biomass and contents of carbon, nitrogen and 
phosphorus . 
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1. INTRODUCTION 

The overall objectives of the LSNDSAT 2 investigation 
No 28 470 "Dynamics and energy flow through, the Baltic 
e'cosystem" was to study the production, transport and 
sedimentation of organic matter- in coastal and offshore 
areas of the Baltic Sea. 

The most specific objective was to develop a technique using 
the synoptic information extracted from LANDSAT 2 registra- 
tions to obtain basic information on ecological parameters 
such as intensity and large scale distribution patterns of 
phytoplankton blooms, particulate matter and effluents from 
urban and industrial areas. 

Such a technique, if successful, would add a new dimension to 
information on the ecology of the Baltic Sea and constitute 
an important tool for the management of this semi-enclosed 
•sea. 

This report will primarily describe an empirical approach of 
how to extract chlorophyll and water quality information from 
two successful registrations in 1975 and how to process distri- 
bution maps of these parameters in the northern Baltic. 

2. TECHNIQUES 

2.1 Test sites and registrations - 

The investigation was carried out as a joint study where several 
Swedish and other Baltic institutes participated. 

The geographical location of the sea truth stations are plotted 
in Fig. la with reference to the respective Institute. A com- 
plete list of participants are presented in Appendix A together 
with sea truth measurements . Fig . 1 b shows the geographical 
location of LANDSAT 2 registrations recived during the investiga- 
tion period of the Baltic sea. 

Table 1 sums registrations, both those of value and those not 
useful due to cloudiness, artifacts and lack of truth data, 
financial or personnel shortage. 
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TABLE 1 Registration statistics over useful and not use- 
ful scenes due to different causes. 


Characteristics 

Number 

Percent 

Cloudiness and/or condensation 
streaks from aircrafts 

21 

35 

Artifacts in the registrations 

7 

12 

Registrations over land 

. 7 

12 

Missing sea truth due to 

6 

10 

geographic location 
other difficulties 

7 

12 

Successful registrations used 

12 

■20 


Total nvunber of registrations 60 

obtained 


•2.2 Sea truth measurements and analysis. 

The sea truth measurements presented in Figure 1 b were taken 
at 10 fixed stations and during cruises in the Stockholm archi- 
pelago and off shore south of Landsort (see Fig. 2),. 

During the off shore cruises in the Landsort area chloro- 
phyll a was continuously measured at 1 m depth .by pumping water 
into the ship laboratory using the in vivo fluorescence technique 
described by Strickland and Parsons (1972) . Figure 3 shows a 
compressed recorder output from the transect August 6. The chloro- 
phyll a values have been calculated from the expression found in 
Figure 4 . 

During the continuous sampling, v/ater was piunped around the 
Turner III instrument and collected in a 100 1 jar correspond- 
ing to 2 naut. miles and subsampled for chlorophyll a, nutrients, 
salinity and phytoplankton counts- For analysis procedure see 
Strickland and Parsons 1972. Secchi disc measurememts were made 
using a 30 cm disc and water telescope. 

This report is concentrated around two registrations, July 20 
and August 6 {Figs . 2 a and 2 b) where sea truth measurements were 
performed synchronous to the satellite in the off shore area, 
while the measurements in the Stockholm archipelago were taken 
3-5 days earlier in July and 2-6 days earlier in August. In- 
vestigations have shown that during this time of the year, the 
water in the different basins have almost the same character 



Figure 5. Compressed recorder output along 
the transect August 6,1975 ( fig. 2b) showrvjg 
the in vivo measurement of chlorophyll a. 



Figure 4- Fluoresconce vs, chlorophyll a during sea truth 
measurements August 1?75* 


for weeks due to low inflow from lake Maelar and calm weather 
(Karlgren and Ljungstrom 1975) v?hich persuaded us to use these 
values in order to expand the range of the sea truth measure- 
ments . 

2.3 Satellite data interpretation and processing. 

As a first step in handling the LANDSAT 2 registrations the 
70 Iran negative products were checked to find out if the scene 
was useful or not. This judgement was built on degree of cloudi- 
ness and quality of the negatives and the presence or synchronous 
sea truth measurements. 

Selected scenes containing relevant information, e.g. algal 
blooms, were further interpreted qualitativly using a light 
table and a stereoscope. 

Direct measurements on the 70 mm products were carried out on 
an automatic instrument, IRIS (Image Reading Instrument System, 
ref. Appendix B) where density differences in the algal accumula- 
tions were recognized. 

This technique has not yet been fully tested but the results 
obtained are promising for a rapid areal quantification of 
specific characters of the surface water, and will be tested 




Chlorophyll a 
measured 


■pigwcs 5* Graf showiiig computed •vs- measured chlcroobyll a 
using stepwise linear regreaGlon* i^ea truth was 
collected off .shore and in StpcJcholsi archipelago* 
LAI?:>3A’:-2 digital numbers in registration August 6, 
1575(13* 2ig6**':C’l720Q) were useC:* 


V = - lS.o6-125j.75xHand 6/5.4-155' .4* banu 6/4 4 ' i;s,8 x bani4/5- 

secchi 



?iguro 6. Graf showing the relation corapute ' vs^ neasured turbidit7 
(secchi disk) using stfcwiae linear regresaion# .^ea tr^th 
was collected off shot? and in Stockholm srchincl-igo. 
liAWDJA'T^.S digital number in regt.^tration August 6, 157; 
(ID. 2196-0917^00) were used. 


further t • 


Quantitative estimations of chlorophyll a an<3 Secchi disc 
measurements have been made using the computerrCompatible 
Tapes (CCT) processed on an IBM 360/75 at the Stockholm 
Computer Center, using programs developed at the Image Process- 
ing Section of the Swedish National Defence Research Institute 
(Appendix C) . 

The following procedures were used in correlating the satellite 
data to sea truth measurements s 

The sea truth stations v?ere located on sea charts, off shore 
by the Decca Navigation System. 

Subscene extraction from the lANDSAT registration over the 
areas of sea truth measurements. 

Symbol coded liner printer sheets covering the subscene 
areas of registrat5.on v/ere used in location of the sampling 
stations. Window extractions, usually 24 pixels in the archi- 
pelago and 240 pixels off shore, were submitted to statisti- 
cal treatment. 

Mean values of the digital numbers in. Band 4, 5 and 6 for 
each sampling station were calculated, using the program 
library. 

2.4. Statistical handling of sea truth and XANDSAT data. 

Stepwise regression analysis were used to find out which band or 
algebraic combination of bands showed the strongest relation to 
sea truth, measurements as dependent variables and digital numbers 
of band 4, 5 and 6 as independent variables»..The calculations were 
carried out on an IBM 360/75 using a Biomedical computer program 
<BMD02E)' . ■ . . 

Expressions to be used in generation of chlorophyll a and tur- 
bidity distribution maps by a PDP 11/40 with an inkjet plotter 
display were also obtained (for exempel see Figs. 5 and 6) 

3. RESULTS 

3.1 Eutrophication monitoring 

XANDSAT 2 registration from August 6, 1975 (ID.2196-0917200) 
over the nor th-v/es tern part of the Baltic, have been used in 
an attempt to generate maps for chlorophyll a and turbidity 
(Secchi disc readings) from the inner part of the Stockholm 
archipelago towards the open Baltic. Stepwise regression ana- 





lysi£5 was nsed. to’ find out which band or algebraic combina- 
tion of bands showed the strongest relation to sea truth 
measurements. The digital numbers of band 4, 5 and 6 were 
used as independent variables and the field measurements . 
constituted dependent variables. By using this technique in 
three steps an equation was established. There was however a 
very small increase in significances compared with using just 
one step. Figure 7 is a representative example. 


STEP MU«a£R X 
VARiAeLH ENTERPO UtoiS 5x6. 

m.TlPl.E fi 0.^360 

S’ .5. ERItOR Of £STo 3S.53.99 

AMAUYSIS of VAFJAKCe 

OF 

ftPCRESSSCN X 

fseSlDUAL 14 


STE? l-fUi'.Br-F 2 
VARIftBie EWEREa 's>eai 6. 

KUtTIftf. R 0.9&00 

STO. EFiRCA OF EST, 32,(>1«5 

AIJAtYSJS OF variance 

DF SUM Of SQUARES 
FfeOscSSICH 2 102546. 137 

SeSfDUAl. la 13H3ii.39l 


STEP HUK3Eft 3 

VARIABLE EKTiREO . 5*4. 

PULTIPLS R 9,97 3 S 

. STJ3. ERF.0R OF EST* 2o.9o44 

AltALVSIS uF VARIANCE 

DF SUH DF SOUARES KEA« SQUARE F 
RFfiRtSSJON 3 1*6310.375 55436.789 66.330 • 

RESIDUAL 12 ■ 10067.230 638,936 

Piffiae 7 . Ha3>^^€tle&4a■fci■ra efeaj>le ohowiTjg sasil i»es»eee in sisengtlx -ot 
rela-tioaship using a ? sta? Xineex xagissAioa eaelysie. 


Table 2 shows theMqhest correlation coefficient (r) and the bands 
used in the statistics. As can be seen, no clear differences- 
appear in the selected bands or band combinations probably be- 
cause of the strong correlation between Secchi disc values and 
chlorophyll a which indicates that the phytoplankton popula- 
tion to a great extent affected the turbidity both in the 
archipelago and off shore during these investigations. 


SOUASE F »AH 
61273.062 76. Ml 

3065.953 


SUK OF SQUARES WEAK SQUARE F ^ATJ J 
154512.000 1545X2.000 98.930 

21865.570 1561.826 
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TABLE 2 Correlation coefficient (r ) and band combination which 
yielded the higheet correlation by stepvjise regression 
calculations * 

Chlorophyll a Secchi disc 



r 

Band 

r. 

Band 

duly 

Archipelago 

.97 

4+5 

.89 

5 

Off shore 

.90 

5+6 

- 

- 

Arch. + off shore 

.95 

4x5 

.81 

4+5 

August 

Archipelago 

.72 

4 

,87 

4+6 

Off shore 

.97 

5+6 

.56 

4x5 

Arch - + of f shore 

.97 

5x6 

.90 

5/6 

July -f August 

Archipelago 

.65 

5+6 

•.76 

5+6 

Off shore 

.96 

5x6 

.46 

4/6 

Arch, + off shore 

.88 

5X6 

.87 

6 


The regression analysis for ail values, digital nmtibers and 
sea truth measurements collected during July 20 (ID, 2179-09232) 
and August 6, show that chlorophyll a is 35est correlated to 
band combination 5X6 and turbidity to band 6. 

The strong significant correlation between the two sets of data 
could probably be enhanced by several parameters in the water 

e.g. yellow substances and other dissolved substances (not 
measured here) , 

Our attempt in generating distribution maps met with diffi- 
culties due to differences in the detectors signals or in the 
calibration. "Work has been done at the Institute for Image 
Processing to overcome the striping by smoohing techniques 
built on the floating mean approach where the output of each 
pixel was ^weighed by the surrounding digital numbers. Figures 
S, and 9 show the uncorrected and corrected products from the 
algal aggregation in Fig. 2fo, 

Another approach was tested in Figure 10 where the combination, 
of bands 5x6 found in Table 2 {August 6, archipelago t off shore), 
were used together with reduction of the scale by taking the ■ 
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FIGURE 8. UNCQRRECTED DCT-DATA (BAND 6) SHOWING BLUE-GREEN 
ALGAL AGGREGATION IN THE SURFACE WATER OF THE 
NORTHERN BALTIC (SEE FIGURE 2B), GENERATED AT 
FOA BY HERTZ INK-8ET PLOTTER OUTPUT CONNECTED TO 
A POP 11/40 COMPUTER. 



FIGURE 9. CORRECTEO CCT-OATA PROOUCT (BAND B) SHOWING THE 
SAME ALGAL AGGREGATION AS IN FIGURE 8, 





FIGURE 10. LANDSA7-2 COMPUTER PRODUCT, SHOWING THE CHLOROPHYLL a OISTRIBUTION FROM THE INNER PART OF 
STOCKHOLM ARCHIPELAGO TOWARDS OFF-SHORE AREAS. THE DIGITAL INFORMATION FROM THE SATELLITE 
HAVE BEEN PROCESSED ON A IBM 75/360 AND CONVERTED TO POP 11/40 FORMAT FOR GENERATION BY 
INK-OET PLOTTER. THE WHITE AREAS ARE CLOUDS, GREEN ARE MAIN LAND AND ISLANDS, AND RED-PINK- 
MAGENTA-BLUE-LIGHT BLUE REPRESENT DIFFERENT CHLOROPHYLL a LEVELS. 

CRED>30 yg/l. PINK 29-10 pg/1, MAGENTA 19-10 jJg/1, BLUE 9-4 We/ 1 and LIGHT BLUE.414 ^g/1) 
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the mean value of 3 line pixels and 4 colxmin pixels as the 
output. The different colours generated on PDP 11/40 con- 
nected to an inkjet plotter, represent different relative 
density classes of chlorophyll a concentrations. 

In Fig. 10 green represents land or skerries, white is clouds 
and red-pink-magenta-blue-lightblue show decreasing values of 
chlorophyll a. 

When a more satisfactory method of over coming the striping by 
separate calibrations of the sensors is available LANDSAT data 
could probably easy be used in productivity monitoring .along 
the coast. 

3.2 Urban and industrial effluents. 

Because of above mentioned difficulties only 70 mm products 
have been interpreted concerning the effluent studies . By 
this crude technique two main systems can be seen, 

I. Effluents having an archipelago as a recipient characterized 
as follows (Stockholms archipelago. Fig. 2) : 

Ko plume or small emmisions can be detected in the LANDSAT 
imagery and there are relatively high ohj-orophyll a contents 
in the water within the area. 

The recipient functions as a sink for nutrients and organic 
matter by sedimentation. 

II. Effluents having no archipelago: (e.g. Wisla-Poland, Venta- 

USSR). . Plumes can be detected outside the river mouth 
(Pigs. 15 and 16) . Nutrients and organic matter are distri- 

A 

buted out into the Baltic. As no truth measurements have 
been reported from these areas no further interpretations 
have been done, but as a subjective impression areas having 
no archipelago could possibly act as nuclec for algal blooms 

3.3 Algal blooms. 

The phenomenon called algal blooms is here understood as accumu 
lation of planktonic algae in the surface water. During the 
summer season the Baltic phytoplankton population shifts from 
a dominance by diatomees and dinoflagellates to biuegreen 
species ( Nodularla spumlgena , Aphanizomenon flosaguae ) and 
monads. See Pig. 11. 
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Figure 11. A mosaic of the 
water below the clouds are 
from August 4-6, 1975 were 


Baltic Sea, where features 
aggregations of blue-green 
used in this composite. 


in the surface 
algea. Registrations 





The aggregation of algae in the surface v;ater is most 
pronounced for Nodularia spumigena . It occurs when the 
algae lose control of their gas vacouls and begin to 
float, ^his is a rapid process, within 24 hours for certain 
species (Reynolds 1971) , but Figure 12 is more relevant for 
Baltic conditions . 
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Figure 1 2 Vertical diatributlon diagrams of Anabaeim cirditalis in Cross Mere, Shropshire, ’ 
- dutitig 1^68. Bo-drawn from Reynolds (1971, fig. 7). 


This pertubation of the buoancy-regulating mechanism is 
caused by several environmental factors such as nutrient 
deficiency, and turbulence decrease. 

In the Baltic the blooms are often observed after a period 
of calm weather, temperatures above 17^C in the surface 
water and high insolation. The aggregation often occurs at 
the borders between Langmuir cells and in the area between 
off shore and coastal waters. 

Figure 13 shows a picture sequence built on LANDSAT 2 regi- 
strations and helicopter observations reference during the 
blue-green bloom in 1975 indicating the irregular appearance 
of the bloom in the northern Baltic . 

This means that total biomass estimations from satellite data 
must be carefully interpreted, although the qualitative in- 
formation obtained has given Baltic ecologists another scale 
to consider in sampling design and measurements from research 
vessels . 





Z*eft: Eesifitration 20,1975* 
documenting the first phase of 
the blne-EJ^een algal assregaiion 
in the surface water* 

Ki^t: Blue-graen algae evenly 
distributed in the siirface water 
(0-2m) starting about 4-*5 naut* 
niles from the nainland* 




Left: Ho algae were detected -in 
the NW^^S route but occurei in 
the ¥-E fligth path. 

Sight 1 SeaTy aggregation in the 
northern pc^t with ch^ging 
deaeitles and directions. The' , 
other route show moderate 
density in the surface water* 



LeftJ 31uo~grecn algae present, in 
the surface water but to a lesser 
ertent* 

Sight: Ii^dlDSAT*^ registration 
August 6,1975 shov;ing the finalo 
Diego of the bloom » 
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Figure 15 « A set of sitaations showi-ctg the variations in the "bl-ae-green 
algal bloom appess^eace in the H\v' Baltic during 1975* Dotted lines indicate 
•helicopter fli^t path* 



3.4 Current patterns 

In Figure 14 th.e LAKfDSAT 2 registration (August 9, 1975 
ID. 2199-09351) shows the main mixing region of saltier North. 
Sea water (ca 30 °Ao) and Baltic water (ca 7-10 *^/oo) . 
Information from the satellite indicates a possible inflow 
situation of water from the North Sea through the Qresund 
where the Baltic water is tagged by the blue— green algae in 
the surface. In this very complex current pattern the arrow 
indicates upwelling due to a coast-parallell winds (E‘ss4 m/sec.) 



Figure 14 LANDSAT 2 registration (10.2199-09351) August 9 , 
1975, showing a part of the main mixing area of 
the North Sea and Baltic waters. The streaks in 
the surface water are' blue-green algae. 


1 % 


Figure 15 is a recording registered on August 6, 1975 
(ID. 2196-03183) , over the southern part of the Baltic sea out- 
side Poland. The structure in the surface water indicates 
heavy . amounts of algae (see arrow) about 6 naut. miles from 
the coast and effluents from rivers point at a SW going 
coastal current. 



Figure 15. LANQSAT-2 registration August 6,1975 ( ID . 2196 -09183 ) 
off-shore 'Poland. Ths configuration in the surface water is 
aggregation of blue-green algae. The arrow indicates ths eff 
luent of the river Wisla. 




Figure 16. LANDSAT-2 'registrations from August 4, 1975 
(ID.2194-09062..and 2194-09055) of the Baltic Sea off-shore 
USSR, showing aggregations of blue-grsen algae in the surface 
water • The outflow of river Vanta can also be noticed. 
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Figure 16 (August 4, 1975 ID. 2194-09062 E 2194-09055) 
again shows a south going coastal current along the USSR 
coast and sharp streaks indicate the borders between the 
watermasses . The accumulations of algae beneath the cirrus 
clouds Indicate a rather calm situation in the central part 
of the Baltic but the upper part indicates north going 
currents. Counter-clockwise gyres can be seen along the 
northern coast of the island of Saarema. 

The last example of surface currents information obtained 
from LANDSAT 2 data is taken from the southern coast of Sweden 
Fig. 17. In this imagery from August 6 (ID. 2196-09181) a- 
nother counter-clockwise gyre is shown. Moreover there are 
indications of upwelling situations along the coast of Sweden 
and off the island of 5land due to K-NE winds. Sea truth 
measurements performed in the area to the lower left showing 
heavy streaks gave chlorophyll a values of 305 ug 1 ^ and 
pheopigment of 399 ug 1 altogether ^700 ug 1 ^ chlorophyll a, 
compared to ca 0.5-1 ug 1 ^ in a water column 0 - 6 m/ which 
indicated a he avy aggregation near the surface. 



Figure 17. LANDSAT 2 registration August 6 , 1975 (ID. 2196- 
09181) . The streaks in the surface water are 
blue-green algae showing a complex current pattern 
on the surface 



There are to date- few estimations of the role of nitrogen 
fixation of Nodularia (Brattberg 1974, Hiibel and Hiibel in 
print) and still fewer calculations of the total amounts of 
biomass and fixed nitrogen in the blooms {Nyqvist 1974, 

SstrSm 1976) . 

The trigger of the blooms are still obscure. They occur during 
the warmest part of the year and it is tempting to point to 
temperature as a major factor. Phosphate as an essential 
element for -the alga could also indicate upwelling areas as 
nuclei for the blooms which during their surface floating 
phase could then be wind-driven (Jansson 1977) . The assumed 
upwelling areas in Figs. 14 and 17 would then act as "nutrient 
windows" letting through bottom generated nutrients to the 
euphotic surface zone . 

Another source of organic materials are the sewage plumes from 
rivers or urban areas. Quantification of the effluents, which 
have been registrated along the coasts of USSR, Poland and DDR 
have not been made because of difficulties in obtaining ground 
truth data. The whole of the Stockholm archipelago is dominated 
by the urban processes which ultimately show up as effects on 
the environment in increased chloropi^fll and turbidity. 

In order to measure the sedimentation rate of the sinking, 
dying algae during the declining phase of the bloom, sediment 
traps had been anchored in the northern part of the investiga- 
tion area. Unfortunately none of these could be recovered, 
probably trapped by trav/ls or fishing nets in these interna- 
tional waters. Although therefore nothing can be said about 
the sedimentation of the registered bloom, previous experience 
has shown that a heavy fall-out from’ Nodularia exists and also 
reaches the bottom (Nyqvist in print) . Due to the sensitivity 
of the deep basins to increased loads of organic matter this 
is of course of great' interest where the decaying Nodularia 

material is deposited. This calls for an intensive survey of 
the blooms throughout their entire phases more extensive than 

the 18 days upon which this investigation is based (Pig. 13) . 
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4.2 The importance of large scale, synoptic quantifications. 
Though water is a fluid medium this by no means secures a 
homogeneous distribution of its suspended solids neither 
living nor dead. Patchiness of plankton populations is a 
well .known phenomenon in the oceans . For modelling of the 
pelagic systems and in order to give a true representation of 
such crucial problems as zooplankton grazing on phytoplankton 
the critical size of the plankton patches has to be detemined 
This has up till now been done almost exclusively by measure- 
ments from ships, at best as continuous tracking with towed 
fluorimeters such as we in some cases employed (Fig. 3) . 
Relying on a small number of fixed stations only - which is 
often the case due to lack of money - results in a high degree 
of uncertainty. The bloom might be concentrated in other parts 
of the sea. 

To date, no information exists on large scale patterns of the 
Baltic plankton which makes a total quantification most un- 
certain. 


4.3. Remote sensing as an operational technique. 

The increasing need for large scale, synoptic surveys in the 
Baltic can only be met with through remote sensing technique. 
Though only a few successful sceneries could be obtained in 
the present investigation these results have already greatly 
enriched the Baltic environmental sciences and induced new 
sampling schemes. The bluegreen blooms appear as an even more 
important Baltic feature which will have to be monitored in 
the future. The remarkably good correlation between concentr- 
ations of chlorophyll a and LANDSAT-2 bands has presented us 
with a powerful tool with which, to survey the archipelagos of 
the Baltic with regard to eutrophication processes. Using 
the Nodularia blooms as "tracer substances", current patterns 
can be studied and critical sizes of gyres> and patches evalu- 
ated. The strait between Sweden and Denmark, the Oresund, also 
also recorded in the LANDSAT- imagery reveals itself as a typi- 
cal mixing area (Fig. 14) , 



4.4, Problems connected with the utilization of LAWDSAT data. 

One of the main drawbacks has been the late delivery of the 
data which in this case severely retarded the whole investi- 
gation. It is not known where in the process of information 
transfer this delay occurred but the first sceneries were to 
hand in July 1975 and in January 1976 the summer registrations 
from August 1975 were received. Also the final date of launching 
of LAHDSAT-2 was for several weeks obscure which made a possible 
survey of the spring bloom of that year impossible. 

Another complication has been the quality of the material. Due 
to either differences in the detector signals or in calibration 
the translation of the bands to digital maps was characterized 
by striping. This had to be compensated for by special smooth- 
ing techniques . 

A great disappointment was also the lack of registration by 
LMTDSAT-2 on the contracted sampling occasion, July 2 1975. 

It was intended to make this passage a primary test period 
incorporating ground truth data at sea with simultaneous 
registration of part of the area by a multi-spectral scanner 
as a part of a cooperative research study by Swedish and 
French space corporations. For some unlcnown reason the track 
of LANDSAT-2 at this period was deflected westward over the 
Swedish mainland and no registration of our test site was 
obtained. - ' 

The last complication has been the difficulty of obtaining 
expected financial support throughout the experiment which 
enforced the running of this project except for the support 
during data processing and sea truth sampling, with only one 
part-time scientist (B.N.) and which has resulted in the late 
presentation of the ‘ quarterly report and much of the data 
still unprocessed. A rough calculation of the total biomass 
of Modularxa blooms is yet to come for example but it is 
hoped that the method of processing the raw data will be even 
faster and easier- 



4.5, Recommendations . 

Although not especially adapted to aquatic use I 4 ANDSAT -2 has 
been of great use in assessing large scale pattterns of import 
ant marine variables . Our recommendations are concerned less 
with, the LANDSAT functioning system than in advice to future 
utilizers. The striping of the tapes was certainly a drawback 
but should be easy to correct. Our main technical experience/ 
however, concerns the handling of the obtained remote sensing 
products. The digital processing of the data, using the tech- 
niques, the reading instr\iment IRIS and programmes developed 
at the Image Processing Section of the Swedish National Defence 
Research Institute (see appendix) can be strongly recommended 
as a simple and fast technique. As stated in sect. 4.4. these 
instruments can produce far more sophisticated results and 
although the ground truth data does not allow more far-reaching 
conclusions, the results as they now stand constitute a con- 
siderable contribution to Baltic science. 

In our opinion no single expression can be calculated for and 
ecological parameter such as chlorophyll which could be applied 
anytime and anywhere in the Baltic. The relation between sat- 
ellite sea truth data has to be determined with, reference 
to the existing phytoplankton population and other natural 
variations .in living components during^ the biological year. 

5. PUBLICATIONS 

The results obtained so far have been published to a very small 
extent in anticipation of. the criticism of this final draft 
report. ERTS-1 products, however have been used in several 
publications, for example Nyqvist 1974, Jansson 1977. In num- 
erous talks and discussion we have used material both, from 
LANDSAT-1 and LANDSAT-2 to exemplify the need for large scale 
synoptic quantifications of biologically active material. 
Examples of such .talks are: 

First sub-regional Eastern-Mediterranean Modelling Workshop, 
Alexandria, 1975. 

Third Soviet-Swedish Symposium on the control of the Baltic 
Sea Pollution, Stockholm, 1975. 



Conunittee Meeting o£ the Baltic Marine Biologists, Gdansk, 
Poland, 1976. 

Joint Soviet “Swedish Expert Meeting, Riga, USSR, 1976. 

6. CONCLUSIONS 

The specific objectives: 1. "To develop a technique for using 
a synoptic view to obtain basic ecological parameters ..." 

2 . "To develop an operational tool for monitoring the Baltic 
Sea ..." have in our opinion been successfully reached. Remote 
sensing technique as demonstrated by LANDS AT registrations 
can be strongly recommended for obtaining large scale patterns 
of phytoplankton blooms, chlorophyll concentrations, discharge 
plumes and surface current patterns in the Baltic area. 

The present results have pointed to the . phytoplankton patches 
as being of much greater scale than previously thought and 
have suggested new areas for both organic discharge and up- 
welling . The future sampling programmes will certainly con- 
centrate on these findings. 
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SEA TRUTH SAMPLING 

Follov;ing institute participated in the experiment. 

Area 1- - 4 , Oresund 

Department of Marine Botany 
Mr Lars Edler 
University of Lund 
Sweden 

Area 5 , Arcona basin 

Institute fiir Meereskunde 
Dr S . Schultz 
WarneiaUnde 
DDR 

Area 6, 7, Hand Bay and E. Gotland 

Fishery Board of Sweden 

Mr 0. Lindahl 

Lysekii 

Sweden 


Area 8, North-western Baltic 
AskG Laboratory 
Mr Bo Nyqvist 
University of Stockholm 
Sweden 

Area 9, Stockholm Archipelago 

Stockholm Water and sewerage Works 

Eng, M. Cronholra 

Stockholm 

Sweden 

Institute of Plant Ecology 
Mr G. EngstrSm 
University of Uppsala 
Sweden 



^0 


Area 10, Aland Sea 

Huso laboratory 
Mr T . Lindholm 
University of Abo 
Finland 

Phytoplankton sampling in algal aggregations along the 
Swedish coast. 

Swedish Meteorological and Hydrological Institute 

Mr .B. Broman 

Norrkoping 

Sweden 

Swedish Coast Guard 

Stockholm 

Sweden 


The following pages contain sea truth data used in this reporu. 
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3.8 


6 in 


6 m 





SEA TRU‘111 DATA, S'LOJJ-j 1( '! >1 AKC'ii U'lU A. y.' 


DATE: 1975-07-17 
S'i’N: 43 

0 in 6 m 12 - m 


Ttiinpeirature 

16.1 

•S ^/oo 

5.34 

©2 ^l^J/l 

9.4 (101) 

PO^-P ug/1 

1 

Tot-P *' 

19 

Ml -N 
4 

NO.^-N 
IMO^-N " 

Tot-N « 

Chi. a ’• 
lot Chi " 
Coliform no/1 
Socchi disk (m) 

CO 

* 


DATE; 1975-07-17 
STKf: 44 


'!’«!.iDerature 
S ^/qo 
O2 Kig/1 
PO^-P ug/1 
Tot-P '* 

NK -N ■' 

]^aD,^-N 
NO--N ” 

Tot,-M ” 

Chi. a " 

■lot Chi ” 
C<.^iii‘omi no/1 
Socciii uj.sk (sO 


0 m 


16.6 

5.39 

9.5 (104) 

1 • 

19 


6 m 


10 ra 


4.1 



SEA TRUTH DATA, STOCKHOLM ARaiU^ELAGO 
DATE: 1975“08“04 
STM: 3 



0 m 

"2 m 

6 in 

Ta>i{)erature 

18.0 

16.7 

11.8 

S ’'^/oo 

2.12 

2.95 

3.95 

O 2 nu/1 

10.6 (117) 

10.2 {110) 

6.5 (63) 

ug/l 

3 


5 

TOt,“P " 

39 


55 


35 


25 

N 02 -N 

10 


12 


165 


415 

Tot.-N " 

700 


800 

Chi. a " 

22,0 

28.3 


Tot Chl." 




Co tiform no/1 

600 

400 

2200 

Secchi disJc &n) 

1.8 




DATE: 1975-08-04 




STN: 4 — 




0 m 

2 m 

6 m 

Ttanperature 

20.6 

19.5 

11.2 

• S °/oo 

3.42 

3.68 

4.56 

O 2 mg/l 

13.5 {157) 

12.8 (146) 


PO^-P ug/l 

3 



Tot-P " • 

41 



•KH,-N " 
4 

260 



NO,-N ” 

4 



NO^-N ” 

1 ■ 



Tot-N *' 

500 



Chl. a " 

28.3 

31.9 


Tot Chl *' 




Colifonu no/1 

100 

200 

1100 

St^cchi disk (m) 

1.1 





SEA TRLTi'H DAl'A, S'lU-'l'HOlM AiClUm^aJC; 
DATE: 1975-08-04 
STK; 5 ■ 


Ttinporature 
S 'Vcx> 

0^ ng/1 
PO^-P ug/1 
Tot.-P " 

4 

N0,-N •• 

NO^-N " 
Tot,-N " 

Chi. d " 

ToL Chi.” 
Colifcaan no/1 
Secchi disk (m) 


0 ro 4 m 8 

21.0 18.2 11.4 

3.84 3.93 4.60 

12,9 (152) 13.0 (151) 2.9 


50 50 50 

1,3 


DATE: 1975-08-04 
STN: 6 



O.m 

2 m 

10' m 

Tanperature ' 

20.8 

19.8 

8.4 

S °/oo 

3.95 

3.95 

5.30 

©2 mg /1 

11.9 (140) 

12.3 (142) 

1.2 

PO^~P ug/1 

3 



Tot-P 

38 



m “N 

ii 

185 



NO^-N ” 

6 



NO.^-N " 

1 



Tot-N ” 

600 



Chi. a ” 

18.6 

16.7 


Chi " 




Colifonn no/1 




Secchi disk (nt) 

2.1 




(28) 


( 11 ) 



so 



SEA TROTH DATA, STOCiOlOm ARaJlPEl.AiJ ) 
DATE: 1975-08-04 
STN: 7 

0 m ■ ^ ^ 

'i’enperature 

21.1 

19.7 

S ^/oo 

4.09 

4.07 

0.y ng/1 

11,9 (140) 

12.2 (140) 

PO^~p ug/1 

5 


•X'Ot.-P ” 

28 



165 


NO 2 .-N 

3 


NO.J-N " 

■it>C.-S " 

1 


550 


Chi. a ” 
Chi." 

i;.o 

12.7 

Coliforra no/1 

50 


Soochi disk (lii) 

3.0 



DATE: 1976-08-04 
SIW: 9 



0 m 

2 m 

TWipecature. 

20.4 

17.4 

S °/oo 

3.80 

3.84 

O 2 ng/1 

12.9 (150) 

12.2 (135) 


po^j-p ug/l 

TOtl-P 

Mi;-N 
<5 

NO^~N 
NO^-N 
Tot-N 
Chi. a 
Tofc Chi 
Colifom no/1 
SGcchi disk (ra) 


_8 

11 - 

4 - 

6 . 

1 

28 
50 
13 
30 0 
450 


10 


8 

4 

5 


ORIGINAL 
OF POOR.QI^^J'^ 


50 

1.5 


200 


600 


m 

7 

,52 

,5 { 64 } 


m 

.1 

.80 

.0 ( 45 ) 



Si 


SEA TRffiH ORTA, aia.'KHOm /UAjim-i-Jl/O.) 
DAIE: 1975-08-04 
SaW: 10 



0 re 

6 m 

10 re 

Teniperature 

20.2 



s ^/co 

9.98 



°2 ” 9/1 

12.3 (143) 



PO^-P ug/1 

4 



at>t.-p “ 

32 




85 



NO^-K " 

4 



NO^-N " 

1 



‘Itot.-N " 

400 



Ghl. a ” 




Tot Chi." 




Golifom no/1 

c 50 



Secchi disk (m) 

1.4 





DATE: 1975-08-04 
SON: 11 

0 re 

2 re 

6 re 

Temperature ' 

20.9 

18.7 

15.1 

S °/oo 

4.06 

4.07 

4.36 

0^ rog/1' 

12.5 (147) 

12.1 (137); 

8.6 

PO^-P ug/1 

4 



Tot-P •' 

31 



NH.-N " 
4 

165 



NO^-N " 

4 



NOj-N •' 

1 



Tt»t~N " 

650 



Chi. a " 

21,3 

21.2 


Tot Chi " 




Colifoam no/1 

50 

50 

50 

Secchi disk (ir*) 

1.5 




(90) 


SEA TIUJTtl DA-EA, STOCKIiULM Al<aiia^i:E AOfJ 
DATE; 1975~08“04 
SaW; 16 



0 m 

2 in 

6 m 

Tesiipecatvire 

20.7 

19.0 

15.5 

S ®/oo 

4.24 

4.24 

4.40 

C>2 

11.0 (129) 

10,8 (123) 

8.4 

iK)^-P ug/1 

4 



I’ot.-P '' 

28 



Nlij-N " 

50 



ko2“N •: 

4 



NO^-N " 

3 

■ 


IVDt.-N ■' 

400 



Chi. -a « 

12.4 

12.8 


Tot Chi.” 




Coiifom no/1 

50 

50 

50 

Secchi disk (rn^ 

2.0 






(89) 



SEA TRUTH JISASUREMENTS SA^^PL3]) 
BY 

I&. G. SSG3TR0M 
IHSfITUTB. of Pr.AHT ECOLOGY 
UHIT3R31TY Of UPPSALA 


SEA TROTH DATA, STOCKHOLM ARCHIPELAGO. 
DATE: 1975-07-16 


Stn. 

Oil. a 

Seech i 

Abs. 480 

A 

19 

2.3 

.143 

B 

23 

2.0 

.174 

C 

28 

1.9 

.221 

D 

27 

1.8 

.233 

E 

29 

1.6 

.297 

P 

28 

1.4 

.310 

G 

27 

1.5 

.343 

H 



.278 

I 

11 

2.7 

.119 

J 

11 

2.7 

.135 

K 

7 

3.4 

.067 

L 

2 

6.0 



M 


i BA TRUTH DATA, STOCKHOLM ARCHIPELAGO. 
LATE: ISTS-Oy-Sl 


Stn* 

Chi. a 

Secchi 

Abs. 480 

A 

8 

2.6 

.070 

B 

22 

1.3 

.140 

C 



_ 

D 

33 

1.3 

.262 

E 

28 

1,6 

.248 

P 

13 

1.8 

.137 

G 

11 

2.0 

.123 

H 

7 

2.5 

.082 

•I 

4 

3.0 

,037 

J 

4 

3.0 

.033 

K 

2 

4.1 

.012 

L 

1 

5.6 

.010 

M 

2 

6.1 

.016 
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SEA TRUTH MEASUREMENTS SAMPLED 
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Sea chaxt over the Uorth-Vfestern Baltic (area 8) showing 
PA.GB track i“or continnos measure of chlorophyll a (doited line) 

sampling stations(fig«xea). July 20,1975* 



SEA TROTH, THE BALTIC SEA, S. LANDSORT, 

DATE . 1375 - 07-20 

Stn, So/oo Tot.Ciil Chi, a Secchi Phytoplankton 


(ug/1) 



1 


to 

1.7 


2 


1.4 

1.1 


3 


1.5 

1.1 


4 


3.0 

2.1 


5 


2.8 

2.0 


6 


1.8 

1.2 


? 


1.9 

1.4 


8 


1.6 

1.0 


9 





10 

6.4 


.8 


11 

6.4 

1.3 

.8 

Stn. 

I 

6.3 

3.3 

.8 

H 

II 

6.4 

2.5 

1.6 

!t 

III 

6.7 

• 2.7 

2.0 

il 

IV 

6.7 

7.5 

6.4 


(a) u^loVl ^ , 

blue-- monads * diatcmees 

Tob, greens ' dinofla^^ 


6.3 

287 

155 

95 

3 

7.3 

449 

301 

126 

32 

7.3 

177 

100 

77 

11 


821 

696 

73 

52 

6 . 5 

1160 

1065 

48 

47 

6.3 

3953 

3777 

102 

74 




?7n 



SEA TRUTH, THE BALTIC SEA, S. LANDSORT. 

DATE: 1975-08-05 

Stn. 3 o/oo tot, chi. chi a Secchi Phytoplankton 

( ug/1) , (m) ' C 10^ / 1) • 

tot, blue- monads diatomeds 
; * greens diaoflag. 


1 

6.47' 

1.3 

.6 

8.8 . 

117 

14 

90 

13 

2 

6.47 

.8 

.4 


143 

38 

100 

5 

3 

6,37 

.8 

.4 


112 

11 

98 

3 

• 4 

6.47 

.7 

.3 

9.2 

61 

14 

46 

1 

5 

. 6.44 

1.3 

.6 


89 

13 

71 

5 

6 

6.47 



9.2 

120 

22 

97 

1 

7 

6.47 

'2.1 

.9 

■ 6.5 

143 

36 

98 

8 

8 

6.33 

.9 

.5 


145 

26 

121 • 

15 

9 

■ 6.18 

.8 

.5 

8.6 

151 

16 

131 

4 

10 

6.04 

1.3 

.6 


178 

23 

154 

2 

11 

6.01 

1.3 

.6 

6.8 

191 

86 

100 

5 

12 

5.86 

1.4 

.8 


359 

92 

252 

1-5 

— 

5.89 

1.4 

.9 

6.0 

307 

44 

250 

13 

14 

5,92 

1.5 

.9 


298 

54 

232 

12 

15 

6.15 

.8 

.6 

9.0 

133 

23 

103 

6 

16 

6.47 




20,6 

77 

124 

5 

17 

6.19 

1.0 

.6 

8.4 

180 

24 

155 

7 

18 

6.4-7 



8.2 

164 

15 

146 

3 

19 

6.47 

1.0 

.5 


142 

4 

136 

2 

20 

6.47 

1.1 

.7 


150 

14 

132 

4 
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Sea chart showing the track and stations 
for Gontinous sampling of chlorophyll a 
6 August, 1975 in north-western Baltic 
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3 o/oo 

SEA TRUTH, THE BALTIC 
DATE: 1975-08-06 

tot.Ghl. Chi. a 

( us/l ) 

f S m 

Seoehi 

(m) 

LANDSORT. 

Phytoplankton 
10^ / 1 

tot, hlue- monads 

greens 

diatoms 4s 
dinoflag. 

1 

6.56 

.4 

.4 

9.3 

95 

16 

76 

3 

2 

6.58 

.9 

.7 

8.8 

276 

5 

263 

8 

3 

6.58 

.4 - 

.4 

9.0 

78' 

7 

70 

1 

4 

6 . 68 

. 6 

.5 

9.0 

154 

19 

129 

6 

5 

■ 6.74 

1.2 

.6 


163 

39 

120 

4 

6 

6.58 

.7 

.5 

9.7 

173 

11 

158 

3 

7 

6.88 

' .5 

.4 

10.5 

107 

14 

87 

5 

8 

6.85 

.6 

.3 


160 

1 

145 

13 

9 

6.88 

.4 

.3 

9.1 

111 

10 

95 

5 

10 

6.75 

.6 

.5 

9.6 

164 

1 

154 

9 

11 

6.75 

.7 

.5 

9.6 

124 

16 

104 

4 

12 

6.40 

.7 

.7 

8.8 

151 

18. 

129 

3 

13 

6.38 

1.4 

.7 

8.8 

145 

8 

132 

5 

14 

6.44 

1.6 

1.2 

6.0 

300 

48 

227 

25 

15 

6.47 

1.2 

.9 

6.0 

229 

50 

. 147 

31 

16 

6.40 

1.6 

1.1 


389 

46 

281 

62 

17 

6.38 

1.9 

1.2 


301 

54 

131 

26 

18 

6.40 

2.3 

1.3 


466 

51 

365 

50 

19 

6.44 

2.3 

1.1 


384 

44 

317 

23 

20 

6.47 

1.7 

1.0 


438 

27 

319 

92 
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IRIS - A n\n AXIS cn>iPARATnR ANTt MTCRnnEN'SITO'lETER USING Tim 
DIFFERENT SCANNING 'jnnES 


N. Aslundj N. von Gersdorff, R. Norherg and 
J.A. Nordin 

Physics TV, Royal Institute of Teclmology, 
100 44 STDCKHOLM 70 


ABSTRACT. IRIS (Image Reading Instrument .SystenO > huilt by SAAB- 
SCANIA in cooperation with Physics TV, KTH, Stockholm, has been 
conceived with the intent of allowing interaction between three 
parties: an operator, a measuring machine and a con^juter. 

Seen from the computer the photographic plate in IRIS is a read- 
only memory, similar to other mass storages. As such it has the 
additional advantage (e.g. compared with magnetic tape) of allo- 
\vijig simultaneous, random access in two dimensions. At the same 
time it is also available for visual observations, by means of 
optical and electronic displays. 

To provide for this twofold role of the plate IRIS is equipped 
with tTO different scanning systems. A small, rectangular region 
around the measuring point is covered- by means of fast scanning, 
with a repetition freq uenc y fitted to the speed of human apper- 
ception. In addition scanning can be performed by moving the 
tables, as is common in most measuring machines. 

The fast scanning allov/s for visual displays on CRT screens of 
the intensity profiles along two perpendicular directions through 
the central measuring point. Iforeover it supplies the measuring 
machine with the information necessary to perform automatic 
settings on objects on the plate, e.g. on stars. Further, it 
makes it possible for the computer to calculate the gradient of 
the density distribution, e.g. in order to track equidensity 
curves. It also offers the capability of fast raster scanning. 

To make certain a high photometric accuracy IRIS has a calibration 
system which relates all measurements to a scale defined by two 
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separate reference levels. Tliis calibration is synchronized witli 
the fast scanning and it compensates for changes in illumination, 
slit width, PM-tube and amplifier gain. 

Accounts are given on the application of IRI.S to different measu- 
ring tasks. 


INTRODUCTION 

An early exaasple of the use of ttMO different scanning systems in 
a measuring machine is given by the so called .Spira3 Reader. 'Ihis 
machine was originally built by prof. L. Alvarez at Berkeley to 
measure bubble chamber plates. Here a fast rotating system pro- 
vides for an efficient use of the measuring time by allowing a 
desired concentration on the objects of Interest. The movement 
from one such object to another is accomplished by moving the 
measuring tables. 

The same basic philosophy has been adopted in IRTS. However, IRIS 
is a general purpose machine, that can be used to measure any 
kind of patterns recorded on photographic plates . Hence the fast 
scanning has been designed to supply very basic information about 
the density distribution around the measuring point, e.g. the 
size and direction of the density gradient. Like tiie Spiral Reader 
this new machine has been designed to work with a computer. 

Further, the idea of making use of interactive procedures has been 
extended and has led to a design, that allows efficient communi- 
cation between three parties: the operator, the measuring machine 
and the computer. 

These objectives sometimes necessitate an unorthodox design, since 
they may lead to conflicting demands. As an example we may take 
the problem of illuminating the plate. To display the plate to the 
operator, a Large part of it has to be illuminated. On the other 
hand the photometric— readings communicated to the computer should 
not be distorted due to sti'aylight. This raises a demand for a 
very small illuminated area. The method chosen in IRIS is to use 
light of different spectral composition for the different purposes. 

Another example comes from the demand that the operator and the 
machine be close to each other, physically, which is in conflict 
with the demand that mechanical or thermal disturbances should 
be avoided. 

In addition there have also been some more general principles 
employed in the design, work: 

1) To acquire the high accuracy (both in positioning and photo- 
metry) rather by automatic calibrations than by large-scale 
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IRIS - A TWO AXIS COMPARATOR AND MICRODENSITOMETER 

use of hi^h precision components. 

2) To make use of the on-line computer as a tool in the process 
of refining the design. 

The design of the machine reflects the cooperative effort of a 
research group and an industry. Thus a basic idea has been to 
create a multi-purpose platfbnn that can be suited to different 
measuring tasks by modular expansions. This platform comprises 
both hardware and software. 

The present exposition serves to exemplify how these different 
objectives and principles materialize in the machine and in the 
use of it. 



Fig.l, The operator controls the measurements from a console that 
includes both optical and electronic displays. The plate is 
mounted into the machine from behind. The outer cover is 
mechanically isolated from the inner parts of the macliine. 


Tm TWO SCANNING SYSTEMS 

The measuring tables, which constitute the slow scanning system, 
are vertical to avoid bending of the plate due to gravitation. 
Actual positions of the tables are measured by separate linear 
digital encoders. 


ORUSINAL PAGE Q 
OF POOR QUALma 
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•IRtS - A TWO AXIS COMPARATOR AND MICRODENSITOMETER 

images pass either of the tv/o fast scanning devices, ivhich 
effectuate X- and Y-scanning respectively. 

Refocussing can be performed during the course of the measurements, 
by moving the detector assemblies. An unintentional, sideivise 
displacement of the slit when doing this \dll introduce no paral- 
laxes, since it will not change the position of the plate image 
in- relation to the image of the index mark. Further the imaging 
is telecentric so that the ma.gnification will not be affected by 
this refocusing. 



Fig. 3. Reference channels for positioning (in the front) and for 
intensity (fiber optics) allow continous coinparisons \dth 
fixed references. The plate is indicated in the figure but 
not the measuring tables. 
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An example o£ the expansion capability of the design is provided 
by the index mark. It can be' replaced by a grid of black squares, 
which define length scales for the fast scajuiing motions in the 
X" and Y directions. This makes possible very fast raster scan- 
ning with a fixed step length. The first grid produced for this 
purpose has a step length of 25 micrometers. 

The -second reference charnel is for intensity. Once every cycle 
of the fast scanning the illumination of the plate is turned off 
by means of a chopper which is synchronised with the movement of 
the vibrating prisms. This defines a reference level for darkness. 
Further, while the illumination of the plate is off, a light pulse 
that bypasses the plate is supplied to the detectors. The magni- 
tude of the resulting electronic pulse relative to the darkness 
level is compared with a fixed reference voltage. The electronic 
amplification is regulated to keep tliis magnitude constant. In 
this way drifts, not only in the lamp but also in the PM tube and 
in the amplifier are compensated for. Consequently the PM tube 
may be driven with a higher current than normally, since slow 
changes of this tube are of no consequences. This benefits the 
d^amic range of the photometric measurements (presently the dyna- 
mic range is just below 4 decades). 

Since the slit is also included in the chain of elements that i? 
encompassed 'by this regulation, the slit width may be changed 
without changing the intensity scale. The optimal slit width can 
thus be set very conveniently by operator interaction, since a 
change of the slit will only change the form and not the size of 
the profile on the CRT screen; 


PROCEDURES TO TEST AND TO REFINE THE MACHINE 

Interactive procedures have been established to verify and refine 
the merits of the machine as a densitometer and as a comparator. 

To test the photometric qualities the operator decides what ob- 
jects shall be investigated, e.g. the fields of a gray scale, and 

directs the machine manually to these objects. The position coor- 
dinates are stored, together with the photometric data from the 

X and Y detectors. These latter values are obtained by integrating 
the transmitted light over rectangular measuring windows. The one 
side of such a rectangle is determined by how long a part of each 
cycle of the fast scanning the detected signal is integrated. (The 
actual .setting is indicated on the CRT screen by a strengthening 
of the profile along that portion) . The other side is determined 
by the width of a preslit. 

The computer mil take over and repeat the measurements an arbi- 
trary number of times. This gives a good opportunity for studying 
the effects of different changes of machine parameters on these 
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measurements r 

This is also done when testing the machine as a comparator. A 
precision grid plate is measured, and the machine makes automatic 
settings on the crossing points of this grid. This is possible 
since the machine has the ability of setting on the median in the 
X- and Y- direction of any intensity distribution on the plate, 
Aslund Ci965)* The (rectangular) areas' of integration employed at 
these settings can be choosen at will. In the present application 
the median setting is used both to make the machine follow the 
mesh lines and to make it perform repeated settings on the cros- 
sing points. 

During these measurements, which allow a separation of the random 
and the systematic errors, the computer also collects data from 
different sensors placed at arbitrary points of the machine, mea-* 
suring voltages 0 / temperatures. A data analysis is perforated to 
reveal any correlation between the systematic errors observed at 
the measurements on the grid and the changes of these parameters. ' 
In this way it has been possible to trace the main error sources 
in the machine and to take appropriate counter measures. Presently 
the emaronnent of the macJiine is not temperature ccmtrolled but 
at night test, with temperature changes vdthin ± 0.5^C, the machine 
has remained stable for several hours within ± 0.5 micrometers. 


APPLICATIONS OF THE MCHINE 

One major application project concerns the measuring and evalua- 
tion of optical spectra. An interactive procedure is employed, 
where the operator decides whether to employ the median setting 
procedure or to scan the line. The latter can be done step by step 
or by continuous sweeping. From the identified reference lines the 
computer calculates a coarse calibration curve, which can be im- 
proved by including more reference lines. Further, from this cali- 
bration curve the wavelength (or wave number) of a line under in- 
vestigation can be displayed immediately , on a data screen. A main 
feature of the procedure is that every determination can be supp- 
lied with a standard deviation, to accompany it during the sub- 
sequent calculations. 

Another application concerns the measuring of spectral reflectances 
of trees on aerial photographs of forests. The individual trees 
are approximated by polygons, defined manually by the operator. 

The computer refinds the same tree on other plates, representing 
recordings in other spectral regions, e.g, to detemnine the rela- 
tive frequencies of different density values within these polygons. 
The results are presented as histograms on a data screen, and the 
aim of the study is to find out what makes it possible for a human 
interpreter to identify different types of trees on anaerial pho- ' 

tograph. The method employed is representative for a basic idea 

or this machine, to mahe possible conversations betiveen operator. 

computer and machine. • ’ 
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INTRODUCTION ' 


: ^ 


Technological advances during recent years have lead to a need 
for research and development of methods in the new field of com- 

» . » ‘ ' 'v ' 

puter aided analysis of pictorial infornnatioTi* An' obvious example 
of such a need for unconventional image processing/photo interpret 
tation techniques is the optimal reconstruction of imagery recor- 
ded by interplanetary space-probes and transmitted in telephoto- 
like fashion to earth- Another example is the attempt to extract 
earth resources and environmental management information from si- 
milar data collected with scanning photometers over a large area 
and digitized in such a manner that each ground areal element is 

recorded in several discrete wavelength bands of the electromag- 
■* * - * ' ‘ * 
netic spectrum* This report constitutes a brief user^s guide of 

a system of computer* programs being developed at the Image Proces- 
sing Section (355) of the Swedish National Defence Research Insti- 
tute (FOA) for the handling, processing and analysis of^ such multi- 
spectral pictorial data and implemented on the IBM 360/75 at the 
Stockholm Computer Center. The fact' that the methods are still 

j 

being developed implies of course that each program description 
is only an account of -the current .status of the system rather 
than a final documentation product- 

A guiding philosophy - so far - has been to work solely in an ext- 
reme form of batch processing, viz* using one specially prepared 
deck of punched cards, a so called source deck, for each separate 
processing task. Naturally the bulk of the’ system subroutines are 
stored on an internal computer memory to simplify the source deck 
handling and to minimize the program initiation and running times. 
Occasionally one has to use a subroutine in punched card version 
- especially for infrequently recurring tasks, recently completed 
subroutines or adhoc subroutines being developed to deal with some 
specific application task problem not already covered by the sys- 


tem* 
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One specific feature - unavoidable for the huge volume, multi- 
spectral pictorial data considered here -- is" the use of an Inter- 
nal label for every image being processed. It is a must for pro- 
gram administration, for simplified data handling and for data- 
type-dependent documentation purposes. Currently satellite (LAND-* 
SAT (4 channels), NOAA/VHRR ,(2 channels)) and aircraft (Swedish 
MSS-75 campaign (10 channels)) data are available at TOA 355 for 
processing. 

The presentation of the program system has been formed as a desc- 
ription, of the specific source deck used for each separate task 
and with the tasks grouped according to their purpose e.g. data 
editing, image generation and so on* In general every separate 
task .source deck description follows a standard layout; 

• short description of what the program does 

• computer listing of the appropriate BLOCK DATA 

• discussion of those variables in BLOCK DATA which may /must 
be altered 

• specification of which subroutines are used 

• line printer documentation log sheet for the actual run exemp- 
lified by the BLOCK DATA. 

Note that the second and last points refer to actual computer lis- 
tings being inserted in the running text and that all discussions 
of these points are limited to the specific task which they exemp- 
lify. 
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DESCRIPTION OF SOURCE DECK 

A typical source deck (figure 2) consists of punched cards with 
the following general structure. 

Job control cards (at several positions in the deck) 

SUBROUTINE CALLS 
BLOCK DATA 
Dataset card(s) 

Comment card(s) (to be printed on line printer) 

The overall program structure is the same in all tasks. A subrou- 
tine PILOT calls three different subroutines; TODAYS, LEGEND and 
CALLS (figure 1). TODAYS prints on line printer the date and LE- 
GEND the text of the user comment cards at the end of the deck of 
cards. CALLS is used to define one specific subroutine for every 
task. 


PILOT 

TODAYS 1 



LEGEND 

1 

CALLS 





task 



Figure 1. The overall program structure. 

In most of the tasks a subroutine called TGHEAD prints the input 
dataset identification head as part of the task documentation on 
line printer. 
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DESCRIPTION OF TASKS 

DATA EDITING 
CCTMERGE 

Merges and reedits the data from a set of four NASA C Computer) 

C Compatible) TCape)s Cfig^re 5) onto a single tape called MERGED 
Cfigure -6)-. AnMnterim earth rotation rectification is performed. 
Also a 256 bytes label will be produced on the single tape. 


BLOCK OATA 

C FOR USE WITH F0A355 TELEGNQSTICS SUBROUTINE LIBRARY 
‘ COHHON/lNCCT/INCCri^) 

REALMS OONAME 

COMRGN/tAQSAM/ KAROO). USn3itLAQSFC(3»6J,0DNAMEl 10) 

C 

OATA INCCr/0*2|3,'i/ 

OATA DDNAME/«CCri^ • t *CCT2<> •,»CCT3^ 

2«MRGFILEl» .•MRGFIL£2^ t •MRGF ILE3* * * MRGF ILE^» . • HRGFILE5* / 
END 


DATA INCCT/0,2,3,4/ INCCTCD-O CCT 1/4 is not available for 

merging. ' * 

INCCTC2)-2 CCT 2/4 exists etc. 


Subroutines in this^ task: 

MERGE CMODIFY, XATICK, INTBYT , BYTINT) 














TASK OOCUHENTATION LOG FOR F0A355 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE (76-07-lA) 

CCT-HERGING, INTERIM EARTH ROTATION RECTIFICATION AND LINE-INTERLEAVING OF THE 

FOUR LANOSAT MSS-BANOS ON A 9-TRACK 1600 BPI MAGNETIC TAPE LABELLED; MERGED 1039-09381 


SC&NE/FkAHE 10 CCT SEQ. DATA RECORD BINARY FRAME ID 
SOOD-IIKMMSON NUMBER LENGTH SAT. CODE 000 HH MM TS 


IOJ9-0UJB1O0 


I 4 


EXPO SURE 

OAifc! ur 


SCENE CENTER 
LAT LONG 


3296 

SATELLITE NADIR 
LAT LONG 


L 39 9 38 1 0 


. BINARY 
STRIP ID 
0 


SOLAR POSITION SAT (REV. 
ELEV. AZIM. OIKINUMB. 

DEGI 


lAT ID 

51182111 

DEFIN. 
PREO. 
ORBIT 


OOlOOlll 


MSS DATA MODE/ MSS ADJUSTED ACTUAL MERGED 
CORRECTION CODE LINE LENGTH RECOROLENGTH 

3240 13640 

MSS DATA ACQUISITION 
DI RECTI ALASKA 
OR (GOLOSTONE 


JIAUGIZ C NVI-21/EOU-00 N N4 /- 2U/E01 1-05 

1039-0938100 3 4 3296 

31AUGI2 C NA7-21/60U-00 N N4 7-20/E01 l-OS 

1039-0936100 4 4 3296 

31AUG72 C N47-21YE0H-OO N N4 7-20/EOl 1-05 

tick hark positions along THE NASA PRODUCED .^‘IMAGERY EDGES (SIGNED INTEGER FRACTIONS) 

TUP EUOE 3705 EOll-OOl -79614 GQI2-OOI -14126 NO48-00I 0 

LEFT LUGE -9446 N048-00= 964 N04'7-30» 11371 N047-00“ 0 

RIGHT EDGE -11970 =£012-30 -7609 =N047-30 2901 N047-00* 13402 N046-30 

9962 EOlO-OOi 3040 :f01,0-i0) -2606 N046-30J 




SUN 



193-0539- 



RECORDI NASAGSFC 



EL45 

AZ143 

-1- 0- 

NASA 

eUTS e-1039-09381- 

R 

N- 


1 V 3;^ 

9 

3a 

1 0 

0 

0 ' 

snaziii 

0 

00100111 3240 



13640 



SUN 

EL45 

AZ143 

193-0539- 

-1- 0- 

NASA 

EfUS E-1039-09381- 

R 

N“ 



9 

3B 

I 0 

0 

0 

51182111 

0 

00100111 3240 



13640 

* 


SUN 

EL45 

AZI43 

193-0539- 

-1- 0- 

NASA 

E«TS E-1039-09381- 

R 

M- 





BorroH EUGe 


0 

0 

0 


HMPf Nnfvnn 


-39oa ieoii-00 -ipa82 uoii-3o 


0 

0 

.6 


it«« ivn ftn lift 




THF ArriiAi^ ronnZn / Ik KECGRUb'D IN ALt FOUR WAVtLENGTH BANOS*. AFTER THE CORIOLI S-EFFECT ADJUSTMENT 

THE ACTUAL GROUND AREA IMAGED IN ALL FOUR k^AVELENGTH- BANOS IS HERE CONTAINED WITHIN: . . , 

I* __ ^ lt3390J TOTAL (LEFT^RIGHD NUMBER OF EQUIVALENT GROUND POINT LOCATIONS IN THE RECTIFYING MARGIN *156 


i 

I / — 

I / 

(2J40, 7) 


— / 
-/ 


i2340»3234l 


APPROXIMATE NUMBER OF LINES WITHOUT MARGIN ADJUSTMENT - IS 


lN^.KChILb2 ?ifFLsT*'L{Nfls‘"' HEKAOECIHAL VALUE Ff. lI.E. DECIMAL 2551 

IN MRGFILE3 The FIRST LINh IS 5«7 AND THE LAST 1172 
IN MKGFILE4 THE FIRST LINE IS 1173 AND THE LAST 1758 
IN MR0FILE5 THE FIRST LINE IS 1759 AND THE LAST 2340 


MERGED SUdSCENE: 
N4/2E110 


LINEJ 

i 


LINES 

2340 


LOCUSl 

1 


LOCI 

3396 


OF HSSBANOS 
ALL FOUR 


MERGED TAPE ONE MUST SPECIFY THE ABSOLUTE COORDINATES BEGINNING AT 1) OF THE LOCATION OF THE 
FIRST GRUUNU POINT ON THE FIRST LINE TO BE TRANSFERRED# THE RANGES FOR LINES AND PIXEL LOCATIONS ALONG THOSE LlNESt AS WELL AS 
THE HAVhLENGTH NUMBERS COMPLETING THE IDENTIFICATION OF THE SELECTED SUBSCENE 

SOMMER-SCHOUL ALPBACH 1976 
PKACriCAL APPLICATION OF REMOTE SENSING 
W0RK:>H0P 4A: DIGITAL PROCtSSiNG TECHNIQUES 
ORhAUG/FUA355 1976 


DATE 76-07-14 
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REPLIC8 

To solve the problem with only two NASA CCTs, each one containing 
two or the four original NASA GCTs. Replicating one of the two 
original NASA CCTs contained on each one of the ^two-set’ tapes* 


BLOCK data 

C FOR USE WITH POA355 TELEGNOSTICS SUBROUTINE LIBRARY 
R£AL*e OF)NAMe, D5NAME 

CCJMMOM/TGLIB/DSNiME, LINEI, LIMBS. LOCUSl •LOCI ♦MSSBI 10) 
COMHON/Li^aSAM/KARD{3) ,LIST(3)tLAC)SFC( 3»6>,0DNAME( 10) 
C 

DATA ODMAME/»ORIGIMAL *t ePLICA 
5ND 


Nothing to change in BLOCK DATA* 


Subroutines in this task: 
Only CALLS. 



O' 


ftfi^-'^ePt.TCrt PftGF 8 

SCEME lOENTl FICil ION LOG FOR F0A355 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE 

SCENE/FRAME 10 CCT SEQ. DATA RECORD lAT 10 BINARY DATA MODE/ MSS ADJUSTED 
SDOO-HriMMSPN NUMBER LENGTH CORRECTION CODE LINE LENGTH 

1039“0'>3sn00 2 A 3296 SI 182 111 39 3290 

2392 RFCDRCS HAVE GEEN REPLICATED 
REPLICATING CCT29 FROM A LANDS iT-COPY-TAPg 


DATE 76-06-09 



17 


SIXPACK 


Converts 7-track NASA LANDSAT tape to 9-track tape. 


BLOCK DATA 

C FOR USE WITH FDA355 TELEGNOSTICS SUBROUTINE LIBRARY 
C 

REAL*8 DDNAME 

COWMON/LAqSArt/KAROO) aiST(3) tLAQSFC( 3i6 1 tODNAME ( iOl 
DATA DDNAHE/* SOURCE ^S»SINK •/ 

END 


Nothing to change in BLOCK DATA. 


Subroutines in this task:* 
SIXPACK(NEGZWO) 





00 


107-SI XPACK PAGE.f r. 16 

\ -s* 

I . -< t 

SCENE/PRAME 1 PENT IF'J C AT t ON LOG FOR F0A3S5 LANDSAT MSS flULK CCT ARCHIVE 
CONVERTING 7-TRACK LANDSAT CCT TO 9-TRACK 


SCENE/FRAME ID CCT SEQ. 
SDOD-HHHMSBN NUMBER 
1039-0931500 I A 


DATA RECORD lAT ID BINARY DATA MODE/ MSS ADJUSTED 
length CORRECTION CODE LINE LENGTH 

3296 SI509701 39 32A0 


EXPOSURE SCENE CENTER SATELLITE NADIR 

date; UT LAT LONG LAT LONG 

31AUG72 C N68-05/E025-54 N N68-03/’E026-00 


SOLAR POSITION SATlREV. DEEIN. 

ELEV. AZIM. OlftlNUMB, PREO. 

OEGI ORBIT 

SUN EL30 AZ167 204-0539- -1- 0- 


MSS DATA ACOUISITICN 
OIRECTlALASKA 
OR IGOLOSTONE 
RECORD INASAGSPC 

NASA ERTS E-1039-09315- R N- 


2342 RECORDS HAVE BEEN TRANSFERRED 

V a JO* KOUSHANEn”^ I9-TRACK» conversion; 1039-09315 NORTHERN FINLAND 


date 76-05-19 
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WIM)OW EXTRACTION 
MERGTAP 

Transfers the 256 bytes label, reedited, and the selected chan- 
nels of a LAND SAT window from an edited 1600 BP I tape to disk or 
tape. 


BLOCK DATA 

C FOR USF WITH F0A355 TELEGNOSTICS SUBRQUTINF LIBRARY 
C 

REAL’*^8 DDNANE,05NA«e 

CaMWON/LAOSAH/KARDBI iLISTI3) *LAQSFC(3.5 )rDONAH£UO) 
COMHON/TGLIB/DSNAMEtL INEI tLINES, LOCUS If LOCI fHSSB( 10) 
COMMON/PLOTTG/LINERSi LI'STEPf LOSTEFf DASH 
C 

data DSNAHF/'STARNBRG*/ 

DATA LINtlfLlNESfLQCUSlfLaCIfHSSB/ 239f 512,1689, $12,4,5,7/ 

DATA LiSTEFfLOSTEP/ l', 1/ 

DATA 0DNANt/*MRGFILEl*,,»HRGFILE2S «MR&FILH3', •riRGFl,E4* , •MRGFILE5* 
2, »DATASETO*/ 

END 


DATA DSNAME/ *STARNBRGV Name of the window transferred. Is ^ 

given by the user. 

DATA LINEl, See the description of the source 

deck, 

DATA LISTEP,L0STEP/1,1/ LISTEP=1 and L0STEP=1 means that 

every- line and every pixel along 
that line (within the selected win- 
dow) is transferred. 


Subroutines in this task: 
MRGTAP (BYTINT jXNTBYT) 


ORIGINAL PAG^ ; 
OE POOR QUAUnf' - 



ro 

o 


12J-MERGTAP PAGE 14 

TASK (JOCUMENl AT I ON LOG FOR F0A3!>5 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE ( 76-07-lS) 
LAN OS AT: 

TRANSFCBRAL OF I039-093ol M5SYCI SUDSCENE FROM MERGED 1600 DPI MAGNETIC TAPE 
EXTRACTION OF SELECTED MSS6ANDS 


SCENE/FR.AMF ID 
SDOO-HHMHSBN 
1039-0938100 


EXPOSURE 
DATE; ur 
31AUG72 C 


SCENE CENTER 
LAf long 
N47-21/EO11-03 


COTS Merged 
0 2 3 4 


MSS DATA MODE/ 
CORRECTION CODE 
OOlOOlll 


ACTUAL MERGED 
RcCDROLfcMGTH 
13640 


INPUT DATASET SUBSCENE: 

LINEl 

LINES 

L0CU51 

LOCI 

NA72tllO 

1 

2340 

1 

3396 

SELECTED SUBSCENE: 

239 

512 

16Q9 

512 

OUTPUT SUB IMAGE: STAKNBRG 

239 

' 512 

1639 

5U 


OF MSSBANOS 
4 5 6 7 

SAMPLED EVERY* 1 LINE AND EVERY I LOCATIONt THUS YIELDING THE 
4 5 0 7 


STARNDPRGtR SEE AND SURROUNDINGS 


SUMMER-SCHOOL ALPBACH 1976 
PRACTICAL APPLICATION OF REMOTE SENSING 
WOR.XSHOP ^a: digital PROCESSING TECHNIQUES 
Ofthf'AUG/FOA355 1916 


DATE 76-07-15 
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CCTAP 

Transfers selected channels of a LANDSAT ■ window 'from a tape with 
data in original NASA CCT format. An interim earth rotation recti- 
fication and a 256 bytes label will be produced. 


BLOCK DATA 

C FOR USE WITH F0A355 TELEGNOSTICS SUBROUTINE LIBRARY 
C 

REAL>^B DUNAME.OSNAHE 

CCHM0N/LAQSAM/KAR0(3) ,Lt ST('3) t LAQSFCI3>6 I .ODNAHE ( 10) 
COMHON/TGUB/DSNAME.LINEl.LINESiLOCCTlTLOCIrHSSBt lOI ' 

C 

DATA LiNEl.UNESfLQCCTi.LOCItHSSB/ 306, 256, 555, 256,6,5,6,7/ 
DATA OONAHE/'DATASETI »,*OATASETO*/ 

END 


DATA LINEl,... 


See the description of the source 
deck. 


Subroutines in this task: 

CCTAP (INTBYT , BYTIMT ,M0DIFY , lATICK) 


OF POOR QGAUE® 



ro 


92<3-CCTAP PAGE Kr - , 

, ’•* 

T^SK UOCUHENTATJQN UOGvfOR- FQA355 MULTRAVEft DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE 1 76-07-13) 
LANOSAT: ^ 


SCENE/FRAME ID CCT SEO.VOATA RECORD BINARY FRAME ID 
SODO-HHNMSBN NUMBER LENGTH SAT. CODE ODD HH HH TS 


1039-0938100 


1 39 9 36 100 


BINARY lAT ID HSS DATA MODE/ MSS ADJUSTED CORIOLISRECT 
STRIP ID CORRECTION CODE LINE LENGTH LINE LENGTH 

0 SI1821II 0, OOlOOm 3240 3396 


EXPOSURE SCENE CENTER SATELLITE NADIR 

date; ur lat long lat long 


31AUG72 C N47-21/60U-00 N N47-20/E0U-05 


SOLAR POSITION SATI REV. 
ELEV. AZIM. DIR|NUH8< 


DEFIN. 

PREO. 

ORBIT 


SJN EL45 AZ143 193-0539 1- 0- 


HSS DATA ACQUISITION 
OIRECTIALASKA 
OR IGOLOSTQNE 
RECOROtNASAGSFC 

NASA ERTS E-1039-09381- R N- 


TICK MARK POSITIONS ALONG THE NASA PRODUCED IMAGERY EDGES (SIGNED INTEGER FRACTIONS) 

TOP EDGE 3705 EOll-OOl -9814 £012-001 -14126 N04U-Q}| 0 m«»hw*««i» q tinnNon.M 0 niiiiitniitf« 

LEFT EDGE -9446 N048-00= 964 N047-3D* 11371 N047-0D” 0 "•"••'••mmii q <ihnhniihu 0 

RIGHT EDGE -U970 =6012-30 -7609 =N047-30 2901 N047-00* 13432 N046-30* 0 «»""»»"« 0 

BOTTOM EDGE 9962 EOlO-OOl 3040 EOlO-301 -2606 N046-30I -3908 lEOtl-OO -10882 lEOU-30 0 «»«•••<"'•■ 


SELECTED UNgl LI§|g LOJ^^J Lg^^ OF M|S8AN0| 

SUMMER-SCHOOL ALPBACH 1976 
PRACTICAL APPLICATION OF REMOTE SENSING 
WORKSHOP 4A; digital PROCESSING TECHNIQUES 
ORHAUC/FOA355 1976 


DATE 76-07-13' 
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CNESTAP 


A window is transferred from a CNES75-tape to disk or tape. Also 
a 256 bytes label will be produced. 


C FDR USE WITH FO&35>3 TELEGNOSTTCS SUBROUTINE LIBRARY 
REAL=^8 OONAMEtDSNAME 

COnHQN/LAQS AM/KAR0(3) ,L TST ( 3 ) 1 1 AQSFC ( 3i6)tDDNAMEE 1F‘) 
CDNMQN/TGLI 6/DSNAME INEl ,L I NES t LOCUS 1 1 LDC 1 1 MSSR f 13) 

C 

DATA OSNAWE/* LARSTAVK*/ 

DATA LI NElt LINE S, LOCUSl »LCCI ♦MSSB/8 75 1 600t38l »L3 2 <4, 5 1 10/ 
DATA DDNAMe/»CNESTAPES < CNESCENE • t • OA TASETO • / 

end 


DATA DSNAME/^LARSTVKNV Name of the window transferred to 

disk or tape. Optional 8 alfanumeri- 
cal characters given by the user. 

DATA LINEl,... See the descirption of the source 

deck. 

Subroutines in this task: 

BYTINT 3 INTBYT 


QJ' POOR 


PA.GB 

qUAITT® 




tN3 

4 >. 


8<1-CNESTAP page 9 

SCENE IDENTIF tCAnON LOG FOR FOAiSG MUtTILAVER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE 

SCENE IDENTIFICATION REG. DATE OT HH MM SS MS 

SU75030A LARSTAVIKEN OA/07/75 12 7 3A 158 

TRANSFERRAL OF A SU8SCENE TO DISK DATASET 

■i ^ 

f 

f % 

^DATE^t6-06-l7 
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WINDOW EXAMINATION 
NUANCES 

NUANCES counts the graylevels of a window and displays the corres- 
ponding histogram(s) . 


C FOR UiE HiTl'i FDA355 TELEGNOSTICb SUBROUTINE LIBRARY 
C 

REAL*8 OONAHE,OSNA«6 

COHHON/TGLl B/DSNAME, L INEl »LIMES , LOCUS liLOLl tMSSBC 10) 
C0MH0N/LA0SAM/KAR013) , LI ST 13) i LAOSFCI 3i6 ) t OONAMEI 10) 
COMMON/ L£VELS/LEV6L<26 o) 

C 

DATA LINEl.LINES.LOCUSitLOCI iMSSB/ 121 6f ZOi l‘«6l »53 . A/ 
DATA CDNAHE/'OATASETI'/ 

END 


DATA LINEl,... See the description of the source 

deck. 


Subroutines in this task: 
TGHIST ( GSALEV {TGHEAD , INTBYT) ) 



6 20— NUAiNCES PAGt 8 

task OOCUMENTATI on log for FDA3!55 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE (76-09-06i 
LAN DSAT; 


SCENE/'FRAME ID EXPOSURE SCENE CENTER 

SDCD-HHMMSCN DATE; UT LAT LONG 

1330-0952300 18JUN73 C N58-41 / E012-23 


ORIGINAL CORIOLISRECTIFIED 
ecu Si i LIWE LENGTH 
1234 12968 


ro 

(J\ 


4^§i ‘•'J'lJ TT'"? 

LEVELCOUNTED SUBSCENE: 1216 20 1461 53 4000 


ACCUMUL. 

PIXEL 

GRAY 

COUNT 

COUNT 

LEVEi 

7 

7 

16 

71 

64 

17 

155 

84 

18 

1 68 

13 

H9 

2Ub 

3 7 

20 

C^3 

^3 

Ail 

55V 

60 

zz 

4:59 

bO 

23 


107 

2-i 

62y 

82 

25 

fZ? 

99 

26 

39 

112 

27 

^01 

62 

28 

y DO 

49 

29 

9^0 

30 

30 

1012 

32 

31 

1022 

10 

32 

1031 

9 

33 

1 

10 

34 

u^r 

6 

33 

1047 

0 

30 

1 uh9 

2 

37 

1J35 

6 

38 

1055 

0 

39 

1050 

1 

40 

1056 

0 

41 

lObii 

2 

42 

1059 

1 

43 

1059 

0 

44 

1060 

1 

45 


X: 


RELATIVE PIXEL COUNT 



DATE 76-09-06 


ORIGINAL PAGE fe 
OF PCX)R quality 
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HISTOGR 

The program selectively counts pixel gray levels within a chosen 
window and displays a histogram. In one of the channels an inter- 
val is specified by the user. Por pixels in this channel with 
gray levels within the ii^terval a histogram is computed* For the 
corresponding pixels in other channels histograms also be 
displayed. Another option is to give a matrix consisting of zeroes 
and onep. The matrix givqs a mask of ones marking the pixels to • 
be included in the frequency count. 


block data 

C FOP USE WITH FQA35S TEL E5N0ST I C S SUBROUTINE LIBRARY 
C 

integer^z mask, select 
INTEGER * 2 THRSHl,THRSH2,MSCNP 
REM<^8 DDNAMEtOSNAME 
COWON/THRSH/THRSHlrTHRSHZtP^SCNR 
COMNCN^VEK/SELECT fHASKt 2500) 

COMMON/ TQU 6/0SWAME,LINEl , LINES t LOCUS It LOCI f MSS B( 10 I 
COMMON/LAQSAN/KARDO } tLTST{ 3) tLAQSFCt 3,G)tODNAME(10) 

COMMON/ LEVELS/LEVEL 12 56 ) 

C 

C 5ELECT=0 (FOR USING THRESHOLDS) 

C SELECT=1 IFOR USINf. »MASKM 

C IN BAND MSCNR IS PIXELS COUNTING ONLY GRAUEVeLS THRSHl ftNP THRSH2 

C 

c 

DATA SELECT/0/ 

data mask/I tO»A*l ,OTA*'l,P,4*l,Ot4*lfOT ItliOtLt 1 / 

DATA line U LINE St LOCUS! rlOCT , MSSB/7 54 , 14 » 2639,1 2 , 4, 5 1 7/ 

DATA THRSHl rrHRSH2,H5 CNR/ 13 t 18,5/ 

C 

C IN M5SBm,tt. ARE THE SAME PIXELS COUNTING AS TH = Y IN RAND HSCNR 
C 

DATA OONAM^/»OAT4SETI V 
END 


COMMON/ VpK/ SELECT, MSK(x) x must not be smaller than LOCI*LINES, 

DATA SELECT/ 0/ SELECTED For thresholds. 

SELECT=1 For mask (punched matrix of 
zeroes s-nd ones)* 

DATA MASK/1, 0,2*1,.. .,1,1/ If SELECT-1, punch 1 if the pixel 

shall be counted otherwise 0. 

A matrix 10110 
1X110 
11110 
11110 
11110 
110 11 

will require MASK of the form indicated. 
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DATA LIHEl,... See the description of the source - . 

deck. 

DATA THRSHl,THRSH2,MSCNR/13,18,Sy THRSH1=13 The lower threshold. 

THRSH2=t8 The higher threshold. 
MSCNR=5 Number of the threshold 
channel. 

T 

Subroutines in this task: 

HISX02(HIST01(TGHEAD, INTBYT) ) 



951-HISTQGR PAGE 8 

TASK DOCUMENTATION LOG FDR F3A355 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE 176-09-151 


LANOSAT: 


scene/frame id exposure scene center 

SODD-HHMMSfiN DATE; UT L4T LONG 


ORIGINAL CORIOLI SRECTIFIED 
CCTISIS LINS LENGTH 


1043-0957400 04SEP72 C N58-33/E011-14 0234 3356 


INPUT DATftSET SUBSCENE; 

LINEl 

LINES 

LOCUSl 

LOCI 

□ F 

MSS0ANOS 

KROPPEFJ 

bOl 

256 

2432 

256 

4 

5 6 7 

LEVELCOUNTED SJBSCEME: 

754 

14 

2639 

12 

0 

5 0 0 


ACCUMUL* 

PIXEL 

GRAY 

COUNT 

COUWT 

LEVEL 

a 

a 

13 

33 

25 

14 

60 

27 

15 

85 

25 

16 

106 

21 

17 

110 

4 

la 


RELATIVE PIXEL COUNT 



lOD PER MILL 

I 

— 227 PER HILL 

— 24S PER MILL 

— 227 PER HILL 
X- 190 PER HILL 

t 

I I 

100 PER CENT 


to 

\D 
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FOA 356 

Causes a window to be .symbolcoded. Every gray level gets a selec- 
ted symbol. The r,esuTt is- displayed on line printer. 


SLOCK DATA 

C FOR USE WITH F0A355 TELEGNQ5TICS SUBROUTINE LIBRARY 
C 

LOGICAL^'l CODE 
REAL*8 DDNArtEfDSNAMG 
COHMaN/COOtycaOE( 256) 

COMMON/LAOSAH/KARDIi) tLI STt 3 ) i LAQSFCt 3 f 6 ) , DDNAM£< LO) 
C0MM0N/rGLIb/D5NAM£vL iNEl tLI NES t LOCUS 1 »L0CI tHSSBl 10) 

C 

DATA LI Ncl, LINES, LOCUSl,LaCt ,MSSB/ 12 1 6,20,iA6i , 53 , A/ 

DATA CODE/15** »t 1* • 1 • , 2 S I* • 3* , !*• A » , 1* »'5S 1* * 6 • , 1* « 7* , * 8* , 

2 ,!*• A* ,1**6*,1**C* rl**0*f !*•£* »1**F*»1**G*, , 

3 I* »I • J* T i**K* , a* r • M* ,*N< ,*0' P*f 'Q* ,*R* f*5*f 

4 ,» V* ,210** '/ 

DATA 00NAME/»0ATASETI •/ * 

END 


DATA LINEl,,,, See the description of the source 

deck. 

DATA CODE/ 15*' ',l*'l',*,. Means that the 15 first gray levels 

, 'u' , 'v' ,210*' '/ are 'printed' with blanks, the 

next with '1' etc. CODE must have' 
256 values. (l*'u' is the same as 

'u')- 

Subroutines in this task: 

TGCODE (ENCODE (TGHEAD , INTBYT) ) 



707-POA356 PAGE 
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TASK DOCUMENT AT ICN LOG FOR F0A3D5 MULTILAYER DIGITAL IMAGt DATA HANDLING AND PROCESSING ARCHIVE (76-09-06) 


L4NDSAT: 


SCENe/f-RAME IC 
SDDD-HHMMSBN 
1330-0952300 


EXPOSURE SCENE -CENTER 
DATE; UT LAT LONG 
18JUN73 C N5B-41/E012-23 


ORIGINAL C ORI OLI SREC T I F I E D 
CCTIS): - LINE LENGTH 
1 2 3. 4 12968 


ll^PUT DATASET SUBSCENcJ 
N58AE122 

GRAYLEVELCODEC SUBSCENE 


LINE! 

LINES 

LOCUSl 

LOCI 

OF 

MSSBANDS 

1201 

128 

1387 

128 

4 

5 6 7 

1216 

20 

^ 146L 

53 

4 

0 0 0 


7A7797677779AA9o9999A?7DA9F“b90EFcA99£HL0HEDAADAAHH00 
8/ia8Ad77b6a7773,7BAA.edbAD0b8b 7,8 DDGGDADK OGEE E6 DDEGE EHDD 
LLb968939b6669CB898G8Bbb9b989b96699aBldlB£EcFECBEINIF 
0 C8 7b b6B6 67 7 67B C 888/ 77 8687 87 db 3436b 7CHJHCBBDFCCCBBJTF 
A fl7AB986A77 7 77A3AB8AA7 77Ab /Aabb53556BCKaLGCCCCCD0GGGa 
9 77CODL)GD97 /779A 90FFF9 79AC V9 964^ 4437ACF JFCC[/0PF6F OGGJ 
y9AaFFFcA967799A9ALV5H99o37644344A69FHPE6EEL)0eFD990EE 
A A0GHGDEHA7 /8AA8AAKSOOK894644A44467 i:GH EDOGHDAEHGGHGDD 
tFLJFIFFF95669B9BBIl EFE5244444444468CECBB8CCCEIECCCCC 
FFFHJHHFCB778B3 V8 88FF7 333-:>33 3 333 3 346BeiBHNFCB88CDCOOOC 
6CGGGCBbbB7AAAAAeA80A73o53233353i5367ACH3HCdCBBbDDO0C 
GG0CCCCDDAA99A9799AC97444 34444433 3 349CCFJJ00CCCCFFF6K 
CAACOOAOAA99AA79D;9D97433A434434344b9DA9AKKOAODcEEEHF 
DaEEAAAAADAe8ADLAA0hA74444444444446o006DaDDt HKHGEOOEE 
CCCCCB99yL98b9G8998CC954444422244246CFEC6bbCECCEFFFII 
CB6CCC8v37BBcl87788b87CDbo3o32 3 33333344bCCCbCCB8CHFCODL 
OCCBCCBrtA77A7bAAAAAAbCDbb55b333333235ADDH0CCCflCCCBDGC 
CAACDAA9y77779AA999AA0DF0 74^4 3333 j 3-»49CC0CCCCCDC0DCCF 
E0HEFEA97777y99A7779tKEEFE0743334443oACADOD£0DDDFKFHH 
AAGEOtA876777d6078a8DM6DOAA8444444446OAAADAAAADDHHGEO 



S^rMBOLCODlNG KSY': 

0 GRAYLEVEL 

1 I l„ 

key: 


15 25 35 45 

-i I I I I I I. 

1234:>6?B9AbCOcFGH[JKLMNaPQRSTUV 


55 65 75 

,_l I ^1 I- 1 


85 95 L05 U5 125 

,.l I- I- I I J— --1 1 I — 




Key: 


O^TE 76-09-06 
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COVCORR 

Computes mean, standard deviation, covariance and correlation. 
The- cor relation coefficients' may aTs’o 'he illustrated graphically. 
One can choose an interval in one channel within which one wants 
to have ‘the statistics quantities. 


BLOCK DATA 

C PUR USt 5JITH F0A355 7ELr.GN0STI CS SUaROUTINF LIriRiRY 

in'.tege:r «2 

nirSCh^^Z THRSHUTHRS^H 2 »MSCNRT 
- OuNAME lOSNA^f ^ 

CO^fM^:WKJMPS/K''^MPS ' 

CC'H^C‘H/LA 0 $A^/KARDi 3 ) ,LIST( 5 J ,L 4 : 5 SPC( 3 %* I , PD^AHE ( 10 ) 
COrtM 3 N/NJAMCf/IM«GEfNUANCEl 2 i> 6 ,n 
CCmHOM/PLOTTG/L IN FftS , L I ST cP, LOSTFP, DASH 
COMMON/ 1 GLIB/DSNAM[:,lINcI ,LI NSS, LOCUS 1 ,LOCI tHSS-JUO) 
CnMHON/THRSH/THRSHl, rH«SH 2 ,MSCNH 
C COflMON/wOisKIN/AREA(KOHPS«^KOMPS + 3 *KOMPS ) 

CC>^H^N/WCBKIN/AK=A(B* 5 ) 

COMMGN/D! SPLA/DISPLA 

c 

DATA KGHF 5 / 3 / 

O&TA THkSHi ,THJ;SH2,M>CNR/13t 20,5/ 

OAfA OS MAMt/ DURBAN */ 

data LI UE 1 1 L I N£S t LQCJ 3 1 1 loci , ^ 55 b/ i) 6 o t bt l^bi ♦ 16 » 5 f 7 / URBAN 

DATA IMAGC/1/ 

DAT a ^AJ^NCE/lt2t3,AT3r6f 7,b, 2*9t2 + o’«^l/ 

DATA LISTPP,LPSTFP/ 6 , I*)/ 

UATA DISPLA/'DO V 
UAT^ D 0 N£iMc/ 'DATASPTI V 
cVJ0“ 


DATA KOMPS/3/ Number of channels. 

DATA THRSH1,THRSH2,MSCNR/15,20,5/ THRSH1=15 The lower threshold. 

THRSH2=20 The higher threshold. 
MSCNR=5 Number of threshold chan-* 
nel. 

Accordinly, the only pixels which 
will be processed in the selected 
channels are those pixels which 


DATA DSNAME/’ URBAN V 
DATA DISPLA/'UO V 


lie within the interval 15 - 20 
in channel 5. 

Name of* the window. To be given 
by the user. 

DISPLA=^DO ’ causes a graphic 
display of the correlation coef- 
ficents. Use DISPIA- ^DONT ’ to 
avoid graphic display. 
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DATA LINEl,. .. 


See the description g£ the source 
deck. 


Subroutines in this task: 

COVCOR (TGHEAD,FLOBYT,LINPRT(MAXIMH)) 


• ORIGINAL PAGE I^j 
c®!poor QUAUTY' 



PAGE 13 

TASK DJCUMcMI ATIGN LOG f-OR FOA35& MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE (76-07-271 
LANOSAT; 

SCEAjr/FRAME ID EXPOSURE SCLNr CENTER ORIGINAL CORI OLI SRtC TI F 1 ED 

SDDD'HHM^SBN DATE? JT LaT LONG CCT(S|; LINE LENGTH 

lw39-OS3/3lO0 3lAUG7.i C N47-2 1 /Eu U-C'? 0 2 3 ^ 3396 

INPUT DATASET SUESCENE.: LINEl LINES L3CUS1 L3CI OF NSSflANDS 

STARNPRG 239 512 1689 512 4507 

SELtCl?D SUB^CENE: iJrt^.-5AL 566 b 1961 16 4 5 07 


::NLY PIXLLS within the INTFRV^L l Ib- MSS&ANO 5) HAYE BEEN USED 


. ^SriHArED Mi:^N and SIANUA^D orvIATUN vectors: 

4 2o374193S+n 7.1721932-/1 

5 Lo 7i22b7'= +:i 1 o5ji<.b9'- + ^J 

7 2.i2d226: + Cl UJ892/ortC9 

ESrUUfrD COVARIANCE (L-lAbCiMAL AND UPPER RIGHT) ANO COWPELATION (LOWER LEFT) HATRICES: 

4 5 7 

4 5c 144.,43“-l 1 5oiYLAK = -tl ?p .^ 9^4 Sc-IU 

5 ‘+ool-t669E~ri 2 ©Ep 9277P+(*' 3o6ood95^ — ri 
7 2o5b4^73:-'l ?• 1 7 b55 L. i 3o523f + :» j 

THE CF Rc PRESENT /iTl VES EXTRACTED FRQvl IHIS MUL TI SPHCTR ALL Y IM6GcD OBJECT FQ^ EVALUATION OF STATISTICS QUANflTIFS IS 124 

7 

4 

4 

4 

4 

4 

4 

5 
b 

b f 

3 ureE*Bee?e?r 

5 CSE'feeeb^Er 

5 C 

7 llllliMli 

7 llllimil 

7 IllillllU 

7 lllllUlll 

7 iiimuii 

7 UlUlilil 

4: 5C'> - 6s*u NamCMETERS 

5: 6ft “ 70t. n:.\OMET? RS 

6: 7CO - auC NANOMETERS 

7J BC-> - UCO hAI^JOMcTERS 


DATE 76-07-27 
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HEXPIXEL 

Prints pixel intensity in hexadecimal foinn of one or several chan- 
nels. Optional versions of subroutine MSSDAT allow lumping in groups 
of all selected channels for each pixel (example 1) or all pixels 
for each selected channel (example 2) respectively. 


Example 1. 


C FOR USE WITH F0A355 TELEGNOSTICS SUBROUTINE LiaRART 
C 

REALMS DONAHEfDSNAHE 

COMMON/TGLIB/OSNAME^LINEltLI NES# LOCUS ItLOCI f MSSB( lOJ 
C0MK0N/LAQSAH/KAR0(31 tLISTO) t UQSFCl J t D0NAM6 U0> 

C 

DATA UNEULINESt L0CUS1«LQCI »HSSB/12^Lt 5f 6t 7/ 

DATA OONAHE/'OATASETIV 

END 


Example 2. 


BLOCK DATA 

C FOR USE WITH FOA335 TELEGNOSTICS SUBROUTINE LIBRARY 
C 

REALMS D0NAME,DSNAME 

COMMON/ TGLIB/DSNAMEtL INEl, LI NESr LOCUS i i LOCI f MS$&( 10) 

COMMON/ LA QSAM/KAR0(3) »LI ST( 3 J ,LAQSFC { 3, S ) , DONAMEI 10) 

C 

DATA LINEl,LlNESiLOCUSi»LOCItMSS6/ 306f 16,1365, 16,4,5,6,7/ 

DATA DDNAME/*DATAScTI V 

END 


DATA LINElj... See the description of the source 

deck. 


Subroutines in this task: 

MS SDAT ( TGHEAD , HEXBYT ) , LAMBDA 


page 12, 

QUArjTK . 



Example 1 


u> 


276-HEXPIXEe PAGE’ 14 

TASK oacUHEHTATION LOG FOR FOA355 HULT ILArER DIGITAL IMAGE DATA HANDLING ANO PROCESSING ARCHIVE C 76-07-121 
LANOSATJ ^ “ 


SCENE/FRAHE 10 ; EXPOSURE SCENE CENTER ORIGINAL 

SDOO-HHMHSBM' '' DATE; UT LAT LONS CCTISl: 

1330-0952300 18JUN73 C N50-AI/E012-23 1234 


CORIOLISRECTIFIEO 
LINE LENGTH 
3240 


INPUT DATASET SUBSCENE; 
ORSJON 

MSSOATA SUBSCENE 


L1NE| LINES 


1201 

1241 


LOCUS 1 
128 138T 

5 1448 


LOCI 

12B 

8 


OF HSS8ANDS 
4 5 6 7 

4 5 6 7 


140B1C0F 15DC1E10 ISDClEll 150C1C10 15OC1C10 140CICOF 140B1AOF 150C1F13 
1400ICOF 140D1C10 160ElC0r 16001CU 14001C0P 140B190F 140ai60F 13001C0F 
140B190F 140C1CU 14081EIU 14051C10 140B1AOE 120A190E I40B190E 140C1E10 
l40C18iQ 14001012 15002012 150C1D10 140C1B0F 140D18CF 140CIBOF 15001D11 
150B1D11 150C1O11 15081FU ISOdlFll 15081D11 150C1F11 ISOBiOll 14081F10 



Example 2. 


401-HEXPIKEL PAGE 14 

TASK OOCUMENTAT IQN LOG FOR FOA355 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE C 76-07-141 
LANOSAT: 



SCENE/PRAME ID 

EXPOSURE 

SCENE CENTER 

ORIGINAL CDRIOLISRECriFIEO 

SD0D-HHMMS8N 

DATE; UT 

LAT LONG 

CCTISJ 5 

LINE LENGTH 


1039-093B10O 

31AUG72 C 

N47-21/E0U- 

00 0 

2 0 0 

3396 


INPUT DATASET SUBSCENE: 

LINEI 

LINES 

LOCUSl 

LOCI 

OF 

N472E11D 


306 

256 

1365 

256 

4 

MSSDATA SUBSCENE 

54 

306 

> 1 < 

16 

1365 

MSS5 

16 

“t >l< — 

4 


HS5S 


>|< \1SS7 


14 1416 1616 L6 141617161 4141610 16 18 000 EO FOE D EO FOEOEOFOE3 0000 EOEOEID 
12131414171713141714171714141418 OBOBOEOEOEOEOBOCOEOEOEOFOEOO 06 12 
1 31 4141714 171 31 7 141417191 B 17 1718 OBOCODOOOEOOOCOCOOOOOOIOOEOOOOOE 
14141314141413141714141819171417 OCOCOCOOO DOEOCOEOOOC030FOFODODOD 
131 31414 1413141414131 314 16 18 181 A OCOBOCOCO63COCOOOOO0OB3OOOOOOE13 
13141313 141 314 14141414161414181A OC OCOBOBO BOBOCODOCOCOCOCOCOEOFOF 
1 61414 16 1414 141 414L61616 13 181617 OEODDB33003EO SOOOOODOOOclOOFOEOF 
17 131314 14 131 21 313 13 141 313 141814 OE OBOCOCOBOBOBOBOAOBOCOEOCOCllOE 
14131 41 41413 14 13 13 L3 13 1417 18 1817 OCOSODOIOCOCOCOCOBOCOBOBOEIOOEOD 
13131413141314141414131413131714 OBOBOCOCOBOBOCOCOBOCOBOBOCOCOOOD 
141314 14131414141414141413141613 OBCBOBOBOCOOODQCOCOCOCOBO BOOOEOB 
1414141314141 31 41414131414 14181 A OCOBOCOCO OOCO BOBOBOCOOOCOCOC 1314 
13131316171313131316161616131317 OCOC080EOEOBOBOCOCOCOEOEOEOCOC13 
13121313131313141414131414131314 OCjAOBOBOCOECBOCOEOCOCOCOEOCOBOC 
18131313141314141717171413141414 OE380COBO CO COODOOEOEOEODOCOCOCOC 
13121313141717171714141714141413 OBOBOCOCOCODOEOOODOCOCODODOOOBOB 


19191A191A1C191A23231C171CIF282A OEOFOFOOOEOFOE101514103E 1014181B 
151 BlCl BlBl C17222 220282 A262 4222 C 0 DOFlllOlO 1 13E15 1 5131816 16 15 151E 
131C1C1E1C1C191C1C1E 2C3A372B1F2C OE 11 121 311OF0F 10101220272318 141E 
iAlElDlE20101AlAlElC2837372A242S 1 112101 2121000101010 1A26251D 1516 
131B101D1S171A131315181F2A202D30 OFOE 101 10D0C3EL 03 EOBOCl 41610 IBLB 
16191516I61818181516191B1B22363E lOOEOCOCODODOEDOOCOOOEll 10162426 
i;i917191A1917171AlrieiC232C2A2C lOOOOEOOOEOFOOOOl 1141313 141A191C 
lC171ClClC1819iBiBlCLB1917222A26 iOOE ill 1121 lOF 1 1 1011110F00141A18 
17171CIC1EIC1CIC1CIA17171E232521 ODOE 12121212111 11 lOFOCOE 12 161714 
lS171ClElCiai01ClAlA IA17171C2224 ODOE 1O121OOE121010OF0EO00C11 1516 
IF1A1BIB14221F101B1BIB1B1B1D1F15 1OOD0EOF0E13111 30FOF IDOFOF 12 lOOB 
18191918101916181819101010102024 1OOEOEOE12OE0D3E0 EOF 1112 1211 1315 
171AI91E1F1A1619U1E2S251E1A191E ODOFOE 1214OE0C000F 12 1316 120E OEll 
1919lB19i:iBlBlE22221E20221CieiC 0 FOE 1 1101210101 2 1 5161212 15 111012 
LE1AL91C1A1C1C1E23231F1E ICIAIEIE lOlOOF 10101 1111 11 6141213 10111212 
IA1710101E22202228201E201EIEIAIC 0 FOCI 212131 51 4161712121211 lOOEll 


SUMMER-SCHOOL ALPBACH 1976) 

PRACTICAL APPLICATION OF REMOTE SENSING 
WORKSHOP 4A: DIGITAL PROCESSING TECHNIQUES 
ORHAUG/FOA355 1976 

DATE 76-07-14 
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PR0JEC2D 

Plots 2-dixnensioml projection(s) of the distribution of the pi- 
xels in n-dimensional space on CALCOMP. (LANDSAT, n=4) 


BLOCK DATA 

C FOR USE WITH FOA355 TELEGMOSTICS SUBROUTINE LIBRARY 
C 

REAL*8 OONAMEiDSNANE 

C0HH0N/LAQSAH/KARD(3) tLI ST(3 ) »LAQSFC( 6 3 , DDNAH6 ( XO ) 
COMMON/TGLIB/OSNAMEtL INEl , LI NES r LOCUS 1* LQC I , MSSB ( 10) 

C 

DATA LINElrLINESrLOCUSlTLXI tMSSB/ 1 1 23i 1421 , 24 7/ 

DATA ODNAME/'DATASETI V 

END 


DATA LINEl,.** 


See the description of the source 
deck* 


Subroutines in this task: 
MFPROJ(TGHEAD) 



fl66-PR0JEC20 P£G£ U' 

TASK DOCUMENTAT tOH LOG POO FOA355 HULTiLAVEIl OfOrTAl IMAGE OAtA IfANOUNG AND PROCESSING ARCHIVE ( 76-09-J6J 


UANOSAT; 


SC£T^.E/fRAMt IQ 

f EXPOSURE 

SCENE center 

original CORTOLISRECUFIED 




SOOn-HHMVSBN 

date; UT 

LAT LONG 

CCTf SJ: 

LINE LENGTH 




1330-0952300 

C 

nsb-a^l/eoiz-, 

23 1 

2 3 4 

3358 




INPUT 0ATAS6T 

suesctwp: 

U NEl 

LINES 

lOCUSl 

LOCI OF 

NSS&ANOS 

KROPPPF J 


1373 

256 

1348 

256 4 

5 

6 

7 

2-D PKOJECreo 

SUeSCENE: 

U49 

23 

1421 

24 4 

5 

0 

7 


LAKE 


KROPPEFJAUt OALSLAMDt CLEAR-CUT INVENTORY tORHAUG.HASTENSONI 


OATE 76-09-16 
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IMAGE GENERATION AND MANIPULATION 
FOA 355 

Generates a window on the line printer with the use of nine prin- 
ter gray levels. 


C FOR*'uil FOA355 TELEGNOSTICS SUBROUTINE LIBRART 

C 

R£AL*B 00NAMe,0SN#ME 

COMMON/LAQSAM/KARDOJ ,LIST(3) ►LAaSFCI3f6 ) lOONANE UO J 
COHMON/TGLI 3/ DSNAHEi L INEl tL INESt LOCUS ItLOC I iHSSB( 10) 
INTEGtR*2 IMAGE, NUANCE 
COMHON/NUANCE/I«ACE,NUANCE( 256,3) 
COHHOH/PLOrTG/LIN£R8,LISTEP,LOSTEP,OASH 
C 

DATA NUANCE/19*9, 2*6, 2*7, 6, 5 , 2*V, 3, 2*2, 2 26*1/ 

DATA LI NEi, LINES, LOCUSl , LOCI ,MS58/i2l6,20, UbI.Sa,*,/ 
DATA LlNERBiLISTEP.LOSTEP/'eACHSlf 1/ 

DATA OD.NAME/ 'DATASET I •/ 

END 


DATA NUANCE/19»9,2*8,2*7,6, ... ,226*1/ The pixels having one of the 


dKiDA? .<./ 

'/s 

DATA LINEl,... 

DATA LINERS , LISTEP , LOSTEP / ’ EACH ’ 


19 lowest digital values (0, 

18) get the printer gray- 
level 9 (corresponding to black), 
the next two get printer grayle- 
vel 8 etc. The last 226 values 
get printer graylevel 1 (white) . 
NUANCE must have 256 values • 

See the description of the source 
deck, 

1,1/ If LINERS has not been assigned 
the value ’EACH^ an interpola- 
tion^ reduction of the number 
of output lines is performed to 
compensate for the rectangular 
form of the line printer charac- 
ters*. . 

If LISTEP=1 and LOSTEP*!, every 
line and every pixel per line 
are printed. 
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Subroutines in this task: 

TGPRNT (LNPRHT (TGHEAD , IKTBYT ,NUANCT) ) 


ORIGINAL PAGE Rj 
OF POOR QUAUTY 



706-FOA355 PAGE 8 ^ 

Si 

TASK 00CUHEN1ATICN LOG FJR FQA35S fULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE (76-09-061 
LANDSATS 


SCENE/FAAME ID 
SOOO-HriMMS8N 


exposure 

DATE; Ur 


SCENE CENTER 
LAT LONG 


ORIGINAL 
CCUSt ! 


CQRIOLI SRECTIFIED 
LINE LENGTH 


1330-0954^300 L8JUN73 C N3d-4i/E0L2- 

23 

12 3 9 

12968 


LINEPRINTER SUBSCENE: 1216 

LINES 

12B 

20 

1461 

LOCI 

128 

53 

OF MSSBAND 
4 5 6 

4 0 0 


f 


EVERT 1 LINESf EVERY 1 LOCI 


DUc TO THE selected LINERATF '«sACH'* THERE WILL 8E 20 LINES AND 53 LOCI IN THE LINEPRINTER IMAGE 

DATA NUANCE/ 19* 9t 2* 8» 2* 7, i* 6» 1* 5» Z* At 1* 3t 

2 2* 2,226* 1/ 


fieBt3ici3iiecsea3«K99»ecfiie(B"fl9i — ees- -leeiae u 

BeBBeacBeBBvsBGcieeetoaitCMieii iei ii ii 

AeAAjBaBBat9AAB9BtCCAA««ii4)Aette8SA A AA- 

AABtAAAACMBteAAAaBBBCtieenCACIIGBeA AAAI-AAAAA - 

aePeAAAAeCIBCe8A6AA6et98DAl3AA«Bllfi«AA AAAAAll A 

• ■■Alll LGBfiiGieai SGiOAB-B^BMIUeSA / All 1 

#€31 esaeCiBBECa 99«IC#HUIH#- ll — l##!— 

6E1 I- Stf#6e3€B SaUUHII##- -U 16- 11 

— ffi###®A4AA aiaiCaiHIMAA-AAAAAAA AAAAA 

-AAfttAACiCC— camMUMinaAAA -AAA§A1A1UA 
1 AAAAAfed336A3AA6iH##UH#ain#9HA AAA AAAI UlA 
lAAAAi leesaeais^&AeaaBuiMiHMuaAA- iaaaa — 

lesi 1 loieosaeesai la# letaaaaaaaasaanasiioee lei i 

i€ — 6e636i6BB6iiGei-tiGuaaaaaaaaHaaii-ici i- -ii— 

AAAAAA#tiiiA#aAi£ateecAA«aaaiaaaaBBBnaA — aaaaa-aa 

AAAAA A#«a; ABaacfCAtGA lAaaaauaaaaniBA AAAAAAAAA -aha 
lAAAAAAAecadaae66eaeAAlABBaaBBBBBaaBBdIl iaaaaaaaai a 

,163 1 166# 3 #66666699066 1 1 - 1 caaBBasnBBB9 AA I A AA AA 1 A1 1 A A- 

i _ I eGaaBeeeaaaaao laaaBBUBMB lei 1 1- u 1 1 — 

'ae -i-eaaBaaaaeaBaaai -u6e9aaBiBuaai666i6B66ii -i 


DATE 76-09-06 
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CALCOMP 


Generates a window on the microfilm plotter with use of the 16 
plot graylevels 10 - 25. 


BLOCK OATA 

C FOR USE WITH FOA355 TCL EG>40STI CS SUBf^DUTIt^c LIBRARY 
R=AL»B DDNA^&t DSNAHE 

COMMDN/LAOSAM/KAROO) ,L I S T< 3 ) , LAqSFC ( 3 , & } , DON AMS < 1 0 ) 
COMMON/rSLIb/OSNAMe, L INgl,LINES,LOCUSlfLO:i ,MSS3( 13) 
CGMMUN/PL0TT3/LINER3, LI STEP, LOST EP, DASH 
INr=GER*2 IMAGciNUAMCE 
COMMON/NUANCc/IMAGFt NUANCE ( 256 , 4 I 
C 

DATA LINEltLINES, LOCUSI tLQCI tHSSB/239 » 51 2t 1609 , 5 1 2t 5 / 
DATA nuance / ll=5>=25f 2^, 22 T 20, 19, 17, 15, 2*1^, 2*1 2, 2 35*13/ 
DATA LISTER, L3STEP, DASH/A, 3, 10./ 

DATA ODNAML/* DATASCTl V 
END 


DATA LINEl,,.. See the description of the source 

deck. 

DATA NUANCE/ 11*25, ... ,235*10/ The 11 lowest values a pixel can 

have get the plot graylevel 25, 
the next get plot graylevel 24 
etc. NUANCE must have 256 levels, 
DATA LI STEP, LO STEP, DASH/ 4, 3, 10. / Each input pixel will he repre^ 

sented by a 4*3 (LISTSP*LOSTEP) 


matrix of output pixels of uni- 
form intensity, 

DASH=10. The length in input pi- 
xels of a line segment and of the 
gap in the dashed line used as a 
frame surrounding the plotted win- 
dow. 


Subroutines in this task; 
TGPLOT (MFPLOT(TGHE AD), FRAME) * 





4^5 


338-CALCOMP PAGE 16 

TASK DD^^'JMENTAriON LOG FOR POA355 MULTIIAVER DIGITAL IMAGE DATA HANOLI'JS ASD PROCESSIvlG ARCHIVE (76-07r23) 
LANOSATs 

center ORIGINAL C3RIOLI SRECT IF lEO 


OF MSS3ANOS 
<►507 
3 5 0 3 


SUMMER-SCHOOL ALPBACH 1976 
PRACTI.CiL application OF REMOTE SENSING 
WORKSHOP AA: digital PROCESSING TECHNIQUES 
ORHAUG/FOA355 1976 






1 I 


1039-^938100 

INPUT D&rASET SUbSCcNEs 
STARNBPS 

microfilmed SUBSCFN?: 

DATA NUANCE/ 1U25* 






U W I k 9 


lime lemgth 


N47-21/EOU-03 0 

2 3 4 

3396 

LINEl LIMES 

LOCUSI 

LOCI 

239 512 

1689 

512 

239 512 

1689 

5U 

l>^:24, 1A22, 1*20, 

2*12,235*10/ 

1*19, 1*17, 

1*15 


date 76-*07-23 
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INKJEXS ' 

Generates a color image on paper or transparency using one or 
more of three possible ink colors. Presently the image is gene- 
rated on tape at the Stockholm Computer Center and subsequently 
imaged at FOA 355* Colors can be assigned as an object identifi- 
cation color code (example 3) or as intensity modulated false 
colors for up to three separate wavelength band images (example 
4);. The three base colors used are yellow, magenta and cyan - 
usjiially in that order for increasing wavelength of the multispec^ 
tral band imaged since the color mixing is of the so called sub- 
tractive type. Resulting color mixtures - in theory - using no 
ink (0) and full (15), ink intensity are as indicated in table 1* 

Table 1. Color generation table for the inkjet plotter. 


^subtractive 

ink 

t 

Y 

M 

C ' 

’theorical^ 
image color 


0 

0 

0 

white 


15 

0 

0 

yellow 


0 

15 

0 

magenta 

intensity 

0 

0 

15 

cyan 


15 

15 

0 

red 


15 

0 

15 

green 


0 

15 

15 

blue 


15 

15 

15 

black 


BLOCK DATA 

C FOR US- WITH F3A555 TELESNOSTICS SUBROUTIWE LIBRARY 
C 

DSNAME, ODNAMt 

COMHON/ T3LI B/OSNAME,LINElr LINES, LOCUS L, LOG I, HSS 3(13) 
COrtMON/LAOSAM/KARDO) ,LIST< 3) ,LAQ5FC( 3, S J , OONAWE ( 1 31 
C0MH0N/PLCirTG/LINER3,LISTEP,L0STE?t0ASH 
INTEGcR^Z IMAGE, NUANCE 
CawMONrNUANCE/rHAGF,NUANCE(256,3) 

COrtMON/MSSr \K/HSSY,HSSC,MSSM 
C 

OATA OSNAME/* I>lKJ6TlMt/ 

DATA MSSY,'iSSCrNS$M/6,3,9/ 

DAfa LI N£1,LINES,LOCJSI,LOCI,MSSB< 1315,133,257,132,3,5,9/ 

DATA NJANC£/ 5^ 15 , 1 A, 13 , 12 , L 1 , 1 0 , 9 , 8, 7 , 6 , 5 , 3 , 2 , 233 ♦ I, 

3 2^13, Ut 13,12,7«^llt3<'10,3*9,2*3,2*7,4+5,3*5,3*if,3, 2, 

4 2 / 2 * 1 / 

DATA LISTEP,LOSTEPOASH/6,6, 13./ 

DATA DDNAMf /OATASETT ' , » DATASET I » , » 0ATA5 cT I ’ , M NKIMAG E* / 
cNO 


ORiGiNAt 'PAGE IS 
OF POOR QUALITY 
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DATA MSSY,MSSC,MSSM/6;3,9/ MSSY=6 Channel 6 is yellow. 

MSSC=3 Channel 3 is cyan. 

MSSM=9 Channel 9 is magenta. 

DATA LIHElj... See the description of the source 

deck. 

DATA NUANCE /5* 15, 14,... See example 4. 

..., 3 , 2 , 222 * 1 / 

DATA LISTEP,LOSTEP,DASH/6,6,10./ Each input pixel will be repre- 
sented by a 6*6 (LISTEP+LOSTEP) 
matrix of output pixels of uni- 
form intensity. 

DASH=10. (has presently no func-, 
tion in this task) . 

Subroutines in this task; 

INKJET (BYTINT , TGHEAD , NUANCT , JETSET ( INTBYT ) ) ' 



388-1 NKJETS PiGE 8 

TASK D3CU'*£NTATI0!< LOG =OR FD&355 

rtSS-75: 


HULTILAYE^ DIGITAL IMAGE DATA HA><DLIM6 AND PROCESSING ARCHIVE 1 76-07-16) 


SCENE IDENTIFICATION PEG. DATE UT HH MM SS MS 

S075333A LAPSTAVltCtN 34/07/73 7 34 159 


INPUT dataset SUBSCENrJ LlNEl LINES 

LARSTVKN 1315 1*3 

SUBSCEN5 SELECTED FOB YELLOW 1315 133 

DATA NUANCE/ 5*l5t l*14t l»13t l*12t 

2 ■ 1* d, 1* 7, 1* 6, 1* 5, 

3 tScI*. 1 / 


LOCUSl LOCI 

257 132 

257 132 

l*llf 1*10» 1* 9t 

1* 4, 1* 3t 1* 2f 


OF MSS9ANDS 
3 3 3 4 5 6 

MSS 6 


3 9 10 - 


INPUT DATASfT SUHStENE: 

UINEl 

LINES 

LOCUSl 

LOCI 

OF HSS3^NDS 

URSTWKN 

131^ 

183 

Zbl 

132 

3 0 3^567 

SOBSCINE SELtCTED f:if CVAN 

13I!> 

183 . 

Ibl 

132 

MSS 3 

DATA HJAHCF/ H*15, 

l*14. 

l*lJt 

1»X2. 

l^lOf 

!♦ 9, 


2 1* «, 
3 231* 1/ 

1* If 

1* ?>t 

I** 5. 

3» 

2# 



input DATASET SUBSCENF: LI N= 1 LIMES 

LARSTVKK 1315 193 

SUBSCcNE SELECTfD F3rl M.IENTA 1315 193 

DATA NUANtS/ 2*l3t l*14t 1*19, 1*12, 

2 It, 2* T, 4* A, 3* S, 

3 222* 1/ , 


L3CJEI L3C1 

257 132 

257 132 

7*11, —3*10, 3* 9, 

3* 4, 1* 3, 1* 2, 


3F NS53AN3S 

333456799 13 
MSS 9 


EACH ORIGINAL »IXEL PLOTTED AS 1 UMIFOBM 6* 6 INKJET PIXEL MATRIX, IHE INKJET MAiE THERfF3RE .HAS 1099 INK-LINES AND 792 INK-LOCI 


DATE 76-1*7-16 


ORiGlNAL page 
OP POOR QUALECa 
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Example 3. 

^5 rKii' ’ 

To do a color image of only^^ppe ^channe^l: of * a classified image 
(wi-th type objects identified by numbers) e*g. a 4 x 4 window: 


0 12 1 
7 10 10 3 

1 1 10 3 

S 5 10 10 


We want the digital values to correspond pairwise to a specific 
color, 0 and 1 (blue), 2 and 3 (green), 4 and 5 (yellow), 6 and 
7 (red), 8 and 9 (light red) and 10 (black). 


According to table 1, blue is a combination of magenta and cyan. 
Using full ink intensities digital level 0 and 1 must have no yel- 
low, but equal amounts of magenta and cyan, (table 2). 

Table 2. 


Digital value 


Ink color 





Y 

M 

c 

Blue 


fo 

0 

15 

15 



0 

15 

15 


I 

15 

15 

0 

Green 

i 


15 

15 

0 

Yellow 


14 

15 

0 

0 

■ 

15 

15 

0 

0 

Red 


16 

15 

0 

15 



15 

0 

15 

Light red . 


18 

[9 

7 

7 

0 

0 

7 

7 

Black 

10 

15 

15 

15 

White 

11-255 

0 

0 

0 


Every column must con- 
sist of 256 values 


Table 2 is then converted color by color to: 

DATA NUANCE/2+0, 6*15, 2*7,15,245*0, MSSY 

4*15 , 6*0 , 15 , 245*0 , MSSM 

2*15,4*0, 2* 15,2*7 ; 15 , 245+0/ MSSC 

MSSY is obtained by assigning 'graylevels^ to column Y and so on. 
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Example’ 4. 

To do an intensity modulated false color image with three sepa- 
rate* channels. 

Do as if simultaneously generating three black and white line 
printer images (task FOA 355) and put them together in NUANCE, 

In NUANCE the channels must be ordered as in the variables MSSY, 
MSSM, MSSC by which one assigns desired channel numbers to the 
appropriate ink color. (Note that the colors in the exemplifying 
BLOCK DATA and in the task documentation sheet are erroneously 
ordered, ) 


SIAL PAGE ^ 
X)R QUAl^rrY. 
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RAMTEK 

Transferral of a window from disk .dataset, to tape for -subsequent 
display on the Ramtek TV-monitor. The Ramtek system allows, gene- 
ration of so called pseudocolor or false color images (color mi- 
xing by the additive principle) . . 


.SLOCK. DATA - . 

C 'FOP use with P0A355 TELESNCSTrCS SUBHOUTtNE LIBRARY 

C . . ■ ' 

RFAL*8 DSfSAME.OONAHE' 

COWMON/TGLI 8/DSNAME,LlNEl ,LI NES» LOCOSl ,L0CI t HSSB ( 10». 

' CONMbN/.UAOSArt/KAPOO) ,LIST(3),LA0SF'C(3t6>.DDNAME( 10) 
COMMON/PLOTTG/LINER8,LI STEP*LCSTEP,DASH . 

•INTEdEPVz I MAGE, NUANCE 
CQr'NpN/^NUANCF/IHA6FtNUANCE(256i3) 
CCMMCN/RAMTEK/LINWOSI 160 J 

c 

CATA LI NEl, LINES, LGCUSl , L OCI , HSSB/ 1390,600, 150 .300 ,5/ 
DATA L!STFP,L0STEP,0ASH/1 ,1, 10./ 

CATA DDNAME/'OATASETI •, 'CATASETO'/ 

END 


DATA LINEl,... See the description of the source 

deck* 

DATA LISTEP,LOSTEP,DASH/1,1,10./ Each input pixel will be repre^ 

sented by a 1*1 (hlSTEP^LOSTEP) 
matrix of output pixels of uni^ 
form intensity, i,e. one input 
pixel gives one output pixel. 
(DASK^IO. has presently no func- 
tion in this task,) 


Subroutines in this task: 

TGRTEK (TGRBW(BYTINT,NUANCT,RAMSET)) 



5G7-RAMTFK PAGE 16 


SCENE 

IDENTIFICATTCN LOG 

FCR 

F0A355 

MS5--75 HULTILAVER DIGITAL 

image 

DATA 

HANDLING 

AND 

PROCESSING 

ARCFIVe 

SCENE 

IDENTIFICATION 


REG. 

DATF UT HH 

MM SS MS 











I 

A 

K 

0 

K E . 










INPUT 

OATASET SU8SCENE: 


LT NEl 

LINES 

LOCUSl 

LOCI 

□F 

MSSBANOS 








SU7501CA 


1390 

600 

ISO 

300 

0 

0 0 0 

5 

Q 

0 

d 

0 

10 

R AMTFK 

B/W SUBSCENGrOATASLTO 

1390 

6C0 

L50 

300 

0 

0 0 0 

5 

0 

0 

0 

Q 

n 


EbCH ORIGINAL PIXEL PLOTTED AS 4 UNIFORM \* I RAMTEK PIXEL NJATRIX* THE RAMTEK IMAGE THEREFORE HAS 6fO LINES AND 300 LOCI 
RAHTFK, CLASSIFIED Jf^&GEt ULTUNA 

DATE 76-0b-25 


€s 

f* s 

u w 


Ln 
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MOVE 


Move a window from a disk dataset to another disk dataset* The 
possibility to halve the values of the pixels exists. 


BLOCK DATA 

C FOR USB WITH F0A355 TELEGNOSTICS SUBROUTINE LIBRARY 

logical*! code 

REAL*3 SELECT 
REAL*8 DONAMEtDSNAME 
COMHON/COOE/CODE{ 256) 

C0HH0N/LAQSAH/KARGC3) fLISTi3J,LAQSFCI3f6it00NAHBll0J 
COHHON/COHPR /SELECT 

COMHON/TGLIB/DSNAHE,UNEl,LINES»LOCUSifLOCI ,HSSB< lOJ 
DATA LI NElt LI NtS, LOCUS! » LOCI *HSSB/ 1150,20, 1250, AO, 5, 7/ 

DATA OONAHE/ ♦DATASET I S *0ATASETOV 
DATA SElECT/'COMPRESS V 

DATA CODE /2* ZOO ,2*Z01, 2*202,2*203 ,2*Z04. 2* Z05 , 2* Z06 ,2* Z07 ,2*Z08 , 
22*Z09,2*20A,2*Z0B,2*Z0C,2*ZOO,2*Z0E,2*Z0F,2*Z10, 2*211,2*212,2*213, 
32*2 14, 2*215, 2*216, 2*217, 2*218,2*219 ,2*21 A,2*Z1B ,2*Z1C ,2* ZiO, 2*ZIE t 
42*21F, 2*220,2*221, 2*2 22, 2*223, 2*224, 2*225, 2*226, 2*227, 2* 228, 2*229, 
52*22A ,2*228, 2* 22C. 2* 220, 2*Z2E,2*Z2F, 2*230, 2*231, 2*232,2*733, 2*234, 
62*235f2*Z36,2*237,2*238,2*Z39f2*Z3A,2*23B,2*Z3C,2*Z3D,2*23E,2*23F, 
72*240,2*241,2*242,2*243,2*244,2*245, 2*246,2*247,2*248,2*249, 2*24A, 
82*24B,2*Z4C,2*24D,2*Z4E,2*Z4F,2*Z50,2*251,2*Z52,2*253t2*Z54,2*255, 
92*256, 2*257, 2*258, 2*259,2*254, 2*256, 2*Z5Ct 2*250, 2*25E,2*Z5F, 2*260, 
12*261,2*262 ,2*263,2*264,2*265,2*266,2*26 7,2*263,2*26 9,2*264,2*268 t 
22*26C, 2*260, 2*26E , 2*26P, 2*270, 2*271 , 2*272,2*273 i 2*274 ,2*275 ,2*276 , 
32*277,2*Z7a,2*Z79,2*27A,2*Z76,2*Z7C,2*Z70,2*Z7E,2*77F/ 

END 


DATA LINEl,.-- See the description of the source 

deck* 

DATA SELECT/ * COMPRESS V If SELECT= * COMPRESS ’ , every value 

will be divided by two. 

Subroutines in this task: 

TGSEL (SELECT (TGHE AD , BYTINT , INTBYT) ) 



7C8-MCVE PAGE 14 

TASK DOCUMENTATION LOG FOR F0A355 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE 176-'09-06I 
LANDSAT: 


SCENE/FRAME IC 
SDDD-HHMHS6N 

1 EXPOSURE 

date; ut 

SCENE CENTER 
LAT LONG 

ORIGINAL CQRIQLI SRECTIFIED 
CCrtSI: LINE LENGTH 


1330-0952300 

18JUN73 C 

N58-41/E012- 

23 1 2 3 4 

3240 




SUBSCENEs 


u15.11 

LOCI 

200 

OF 

4 

H|SBAND^ 

TRANSFERED SUBSCENE 

llbO 

20 1250 

40 

0 

5 0 7 


ALL VALUES OF THE PIXELS HAVE BEEN DIVIDED BY 2 
DATE 7b-09“06 


Ln 
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COMPOUND 

Calculate the difference or the sum between two channels. For 
other expressions one must insert the appropriate program state~ 
ments in subroutine OPTION. 


C 

C 


c 


c 


BLOCK DftTA 

FOR use WITH F0A355 T6LECN0STICS SUBROUTINE LIBRARY 

LOGICAL*! OIVTAB,DIYGO 
LOGICAL*! B1W0RK(512) 

REAL*8 DONAHEfDSNAHe 

COMMOWTGLie/DSNAME,! INElf LINES* LOCUSlf LOCI tHSSBt 10 J 

COMHON/LiQSAM/KARDOJ .L1STI3) »LA0SPC(3*6l rODNAHEUOI 

INTEGER*2 IH!fIH2fSIGM 

COHNON /EXPR/ IM1»SIGN,IW2 

COMMON /tab/ 0IVTAB(414>,DIVG0t98J 

C0MMGN/W0PK/81W0RK 

DATA LINE!, LINES, LOCUS!. LOCI, HSSB/1L03, 61. 1291, 100. 5 *7/ 

DATA IMl,SICN,IM2/7,’-',5/ 

DATA ODNAHF/ • DATASETI » , • OATASETO */ 

DATA DIVTAS/2*Z0O,2‘'ZO!,2*ZO2, 2* Z03,2*ZO^,2*Z05,2*2J6, 2*237,2*208, 
22*ZO9,2*Z04,2*ZOB,2*Z0C,2*ZOD, 2*ZOE,2*ZOF,Z*Z10,2*Z!!-,2*Z1Z,2*Z!3, 
32*Z14,2*Z!5,2*Z16,2*Z!7,2*Z18,2*Z19,2*Z1A,2*Z13,2*Z!C,2*Z10,2*Z!E, 
42*Z!F,2*Z20,2*Z21,2*Z 22, 2*Z23, 2*Z24, 2*Z25,2*Z26,2*Z27 ,2*Z23 ,2* Z29 , 
52*Z2A,2*Z2B,2+Z2C,2*Z2D,2*Z2E,2*Z2F, 2*230,2*231, 2*232,2*233, 2*2 34, 
62*Z35,2*Z36,2*Z37,2*Z38,2*Z39,2* Z3A,2*23B,2*Z3C,2*Z3D,2*Z3E,2*Z3F, 
72+240, 2*241, 2* 242, 2*243, 2*244, 2*245, 2*246,2*247, 2*248, 2*249, 2* 2 44 , 
82*Z48,2*Z4C ,2* Z4D,2*Z4E,2*Z4F,2*Z50,2*Z51,2*Z52,2*Z5 3. 2*254,2*255 , 
92*256, 2*Z 57, 2*Z58»2*Z59,2*Z5A, 2+Z5B ,2*Z5C,2*Z5D,2*25E ,2*Z5F,2* 260 , 
12*26 1, 2*262, 2*263, 2* 264, 2*265, 2*266, 2*267, 2*268, 2*259, 2* 26 A, 2*Z6B, 
22+Z6C, 2*260, 2*Z6E, 2*2 6F, 2*270,2*271, 2*272, 2*273, 2*274, 2*275, 2*275 , 
32*277, 2*278, 2* 279, 2*Z7A, 2*27B,2*27C, 2*27 D,2*Z7E,2*Z7F, 2*280, 2*2 81, 
42*262,2*233,2*284,2*285,2*266,2*287,2*288,2*289,2*23 A,2*Z8B,2*28C » 
52*28 D, 2*Z8F, 2*Z8F, 2*290, 2*291,2*292,2*293, 2*294,2*295,2*296,2*297, 
62*Z98,2*Z99,2*Z9A,2*29B,2*Z9C,2*290,2+Z9E,2*Z9F,2*ZA0,2*ZA1,2*2A2, 
72*ZA3,2*ZA4, 2*ZA5,2*ZA6,2*2A 7, 2*ZA8,2*ZA9,2*ZAA,2*ZA8,2*ZAC,2*2AD , 
82*2AE,2*ZAF, 2*280 ,2* 281,2*282, 2*283,2*20 4, 2*265,2*286,2*287,2*288, 
92*289, 2 *ZBA, 2*zaB,2*2BC, 2*280, 2*-2BB,2*ZBF,2*2C0,2*ZC I, 2* ZC2,2*ZC3, 
12*2C4,2*ZC5,2*ZC6,2*2C7,2*2C8,2*ZC9,2*ZCA,2*ZCB,2*Z:C,2*ZC0,2*ZCE/ 
DATA OIVGO/ 

22*ZCF, 2*200, 2*201, 2*202, 2*203, 2*204, 2*ZD5,2*ZD6, 2*237,2*208,2*209, 
32*20A,2*2OB,2*ZDC,2*2O0,2*ZDE,2*Z0F,2*2EO,2*ZEl,2*ZE2,2*2E3,2*2E4, 
42*ZE5,2*ZE6,2*ZE7,2*2E8,2*2E9,2*2EA,2*ZEB,2*ZEC,2*ZE0,2*ZEE,2*ZEF, 
52*2F0,2*ZF1,2*2F2,2*ZF3,2*ZF4,2*ZF5,2*ZF6,2*ZF7,2*ZF8,2*ZF9,2*2FA, 
6Z*ZF8,2*2FC,2*ZFD,2*2FE,2*ZFF/ 

END 


See the description of the source 
deck. 

IMl and IM2 mnnber of channels. 
SIGK='-' or 

/7,'-',5/ means therefore channel 
7 minus channel 5. 

origmalpageB 

OF POOR 


DATA LINEl,.. . 

DATA IMl,SIGN,IM2/7, ’-’,5/ 



Subroutines in this task: 


COMPOZ(TGHEAD,BYTTNT, OPTION) INTBYT) ) 



Ln 

ON 
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TASK DOCUHENTATION LOG FOR FOA3S5 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE 176-09-02) 
LANDSAT: 


SCENE/FRANE 10 EXPOSURE 
SDDO-HHHMSaN DATE; UT 

1330-0952300 18JUN73 C 

INPUT dataset SUBSCENE: 
ORSJdN 

MSSDATA SUBSCENE 


SCENE CENTER ORIGINAL CORI OLI SRECTIFIED 
LAT LONG CCTIS); LINE LENGTH 


M58-A1/E012-23 I 

LINEl LINES 

1073 128 

1103 61 


3 4 

3240 

LOCUSl 

LOCI 

1240 

200 

1291 

100 


OF MSSBANDS 
4 5 6 7 


DATE 76-09-02 
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TRAINING AND CLASSIFICATION 


TYPEHEAD 

Initiating a type-object library. 


BLOCK DATA 

C POR USF WITH FOA355 TFLEGNOSTICS SUBROUTINE LIBRARY 

c ■ 

REAL+8 ODNAMEfOSNAHE 

COMMON/ TGLIB/DSNA«^,LlNElf LINES, LOCUS I » LOCI tMSSB[ 10) 
COMMDN/LAQSAM/KARD'O) rLIST(3) , LAOSFC ( 3,61 , DDNAME tl 0) 
COMMON/ KOMPS/KOMPS 

COMMON/MSTyP=/ MSTYPE ,MSCHIN,MSCMAX 
C 

DATA DSNAMF/’ N584E122 »/ 

data MSTYPE,MSCMIN,MSCMAX/UAr7/ 

DATA KOMPS /3/ 

DATA LT NE 1, LINES, LOCUS I , LOCI ,MSSB/ 1073,2 56, 13^3, 2 56, ^,5 , 7/ 
DATA DDNAHE/ 'DATASETO*/ 

END 


DATA DSNAME/'N584E122V DSNAME is the name of the scene 

to be classified (constructed in 
a prescribed manner from the geo- 
graphic coordinates - here N58- 
4L/E012-23 - for the scene center). 
DATA MSTYPE,MSCMIN,MSCMAX/1,4,7/ MSTYPE-1 For LANDS AT* 

MSTYPE=2 For MSS-75. 

MSTYPE=3 For VHRR. 

MSCMIN Loi/est channel number. 
MSCMAX Highest channel number. 
DATA KOMPS/3/ -KOMPS Number of channels. 

DATA LINEl,... See the description of the source 

deck. 


Subroutines in this task: 
BYTINT 




Ln 

00 


9 84'Tlp,YPEHtAO PAGE 

TAS'k-’^DOCUMENTATION log for F0A355 MOLTILAYER OIGtTAl IMAGE OATA HANDLING AND PROCESSING ARCHIVE (76-09-L4J 


INITIATING TYPE-CLASS TRAI NING-LI BRARV FILEsN584El22 
KRQPPErJALL, DAL5LAN0, CLCAR-CliT INVENTORY lORHAUG.HASrENSON) 


OATE 76-09-14 
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TYPELIB 

Stores original digital data from type-objects for further use in 
maximum likelihood classification* The type-object does not have 
to he rectangular, not even connected, but can be any number of 
arbitrarily formed areas chosen by using mask(s) or produced by 
using thresholds in one channel as determined with the aid of the 
HISTOGR task. ’ ‘ 


C 

C 

C 


r 

c 

c 

c 

c 

c 

c 

r 

c 

c 


BLOCK HAT A 

FQR USE WITH FOA355 TELEGN'GSTICS SUBROUTING LlBRARy 
LOGICAL*] PIXEL 

IN PIXEL(X), X MUST NOT 6E LESS THAN LI NES*LOCI * (NUMBER OF CHANNELS) 
COMMOM/ PIXEL/ PIXEL(2500) 

INTEGEft^a ODJNRfARtANRrCODE 
CONHON/PARAP./ OB J NHL , CD JNRt ARE ANRy CODE 
INTEGER*? SELECT, MASK 
COMMON /VEK/ SELECT, HAS K { AOO ) 
iM^EGeR*? THRSHltTHRSHZ,MSCNR 
COMMON /THkSH/ TMRSHl ,THRSIP ,MSCNR 
REAL *8 DDNAI1E,OSNAHE,OBJNHE 

COMMON/ TCLI e/iliSMAMe,L rwHl ,U NFS, LOCUSl,LOCT ,MSS8l in» 
CQHHON/LAQS&M/rARO{ J) ,L I ST ( 3 > , LAOSFC ( 3 16 ) ,OONAME( 10 ) 

SELECT=0 <F0R USING thresholds ) 

SELECT=1 (PCJR US^NC MASK) 

OF THE TyPE-OBJECT TO (5F ENTERED 
AREANR-nUNBER OP THE TYPE -OBJECT 5U3AKLA TO BE ENTERFD 
C0nE = l (NEW TYPE-CIASS ANO NEW AREA) 

CODE = 2 (NEV; AREA) 

CODE-3 (DCl VZ AN ARCA ) 

C0i>E"^ (UElETE A TYPE-CLASS) 

data select /o/ 

0\ TA MA SK/5* 1,0, I ,0/ 

DATA ODJNflE,OBJNR, AREANR, CODE/ 'P INF • *5,2, ?/ 

DATA LINE! , I TNES, lOCUSl finer ,M5SB/ 1297, ^,1453, 4,4,5,?/ 

DATA THRSHJ ,THRSH2,MSCNP/J, 2*>5,4/ 

data DDNAHS / » LIRRARY * *TEMPQRAPS ' dataset I* , » PR intent V 
END 


COMMON/ VEK/ SELECT, MSK(x) x must not be less than LINES* 

*L0CI when SELECT=1. 

COMMON /PTXEL/PIXEL(x) x must not be less than LINES* 

*L0CT* (number of channels) . 

DATA SELECT/0/ SELECT=0 Por thresholds in channel 

MSCNR. 

SELECT=1 Por mask (punched matrix 
of zeroes and ones). 


ORIGINAL PAGE IS 
OF POOR QUAUnS 
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DATA MASK/5* 1,0, 1,0/ If SELECT=1, punch 1 if the pixel 

shall be counted, otherwise 0. 

A matrix 1111 
10 10 

corresponds to ^DATA MASK of the 
form 5+1, 0,1,0* 

DATA 0BJNME,0BJNR,AREANR,C0DE/ 

/’PINE ’,5,2,2/ OBJNME Name of the type-object to 

be entered- 

OBJNR Number of the type-object 
to be entered. 

AREANE. Number of the type-object 
subarea to be entered- 
C0DE=1 New type-object and new 
. area. 

C0DE=2 New subarea. 

C0DE=3 Delete an area. 

C0DE=4 - Delete a type-object. 

DATA LINEl,./’. See the description of the source 

deck. 

DATA THRSHl,THRSH2,MSCNR/0,255,4/ THRSH1=0 The lower threshold. 

THRSH2=255 The higher threshold. 
MSCNR=4 Number of the channel. 
With these threshold parameters 
all pixels within the selected win- 
dow are accepted as belonging to 
the type object. 


Subroutines in this task: 

REORG (INTBYT , STATIN (INTBYT , TGHEAD) ) 



ORIGINAL PAGE IS .. 
OP POOR QUALITIS 
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TASK DQCUMFNTiT TOM LOG FDR FOi355 MULTILAYER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE (76-09-14) 
LANDSAT: 

SCFMF/FRA^’F ID EXPOSURE SCENE CENTER ORIGINAL CDRIOLI SRECT! FIEO 
SDDD-HHMMSBM DATE; UT LAT LONG CCT(S): LINE LENSTH 


1330-0952300 

18JUN73 r 

M58-4I/P012-23 

1 2 3 

4 

3358 







INPUT DATASET SUBSC^NF: 

LIMfcl 

LIMES 

LOCUS! 

mci 

OF MSSBANOS 





KROPP*^F J 


1073 

256 


1348 

256 

4 5 6 

7 





TYPE LIBRARY: N584F12? 






4 5 0 

7 

□ BJS: 5 

KOMPS; 3 


QBJFCT 

PIXELS 

PBJFCT 

^ PC AS 

LlNFl " 

L INFS 

LOCUS! 

LDCI 

THRES- 

THR 

= S- 

THRESHOLD 

PI XEL S 

AREA 

NAME 

^^'KOMPS 

mumper 






HOLD 1 

HOLD 2 

BAND 

H 

NUMBER 

CLPARC73 

12^ 

] 

1 

1301 

9 

1467 

8 

22 


25 

4 

129 

1 

CLEARr?2 

66 

2 

1 

1292 

4 

1445 

8 

i * 22 


24 

4 

66 

1 

REGROWll 

69 

3 

1 

1303 

6 

1400 

5 

> 22 


26 

7 

69 

1 

RFGP0W16 

60 

4 

1 

1302 

5 

1475 

5 

I 18 


23 

7 

60 

1 

PINT 

96 

5 

2 

1291 

A 

148 3 

A 

0 


2 55 

4 

48 

1 





1297 

4 

1453 

A 

r 0 


255 

4 

48 

2 


KQQPPEFJALLf OALSLANDr CLE6R-CUT INVENTORY ( OP HAUG » WASTE NS ON) 
OATF 76-09-14 


o^ 
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KEYHEAD 


Initiating a type-object key-library. 


BLOCK DATA 

C FOR U5^ WITH FOA35‘? leLECNCSTICS SUBROUTINE LIBRARY 
C 

REAl'«^8 ODNAMEf OSNAHE 

COMHON^TGU D/OSMAMEtLINEI, lines, LOCUS 1, LOCI pMSSBI 10) 
C0MM0N/LAQSAM/KARP(3)fLIST(3J#LAQSFC( 3^6 ) ,OONAHE 110) 
COMMOM/KOHP S /MS TY PE t KOMPS 
C 

DATA DSNAME/»N5^)AE122V 
DATA MSTYPE ,KDMPS/l,3/ 

DATA MSSU/"t,5»7/ 

DATA C0riAM£/«KEYL IB */ 

EMO 


DATA DSNAME/»W584E122 V 
DATA MSTYPE, KOMPS/ 1,3/ 


DATA MSSB/4,5,7/ 


DSNAME is the name of the scene 
to be classified (cf. TYPEHEAD). 
MSTYPE=1 For tANDSAT. 

MSTYPE=2 For MSS-75. 

MSTYPE-3 For VHRR. 

KOMPS Number of channels. 
MSSB(l),... Selected channel numbers. 


Subroutines in this task: 
BYTINT 



2A6-I^EVHEftO PACE- 8 

task DOCIIMEMTATI i:n tor, fop. F0A3';5 multilayer OIGITA-I. image data handling and processing archive ( 76-09-lAI 

INITIATING TYPE-CLASS KEY-LIBRARY F I LE :N5B AE12 2 
KfiOPPEFJALLt DALSLANOi CLEAR-CUT INVENTORY ( ORHAUGtWASTENSONI 
DATE 76-09-lA 


O' 

Lo 
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INTRAKEY 


Evaluates initial type-object statistics quantities in a library 
called KEYLIB to be used for maximum likelihood classification. 


BLOCK OAT& 

C FOR USE WITH FO&355 TELEGNOSTICS SUBROUTINE LIBRARr 

INTEGER*A OBJEC 
REAL*8 ODMSMc,EPS,OBJNME 
REAL*e DELETE 

CONMON/CLASS/ NR0BJ,0BJEC(24) 

COMMON/ DELETE/OELETEt 24) 

COMMON/ESCALA/EPStlTFRB 

f 

C COM«nN/KFYLIB/KEYSfL=NGTH) 

C LENGTH MUS"^ NOT BE LFSS THAN THE BLOCKSIZE OF THE CATALOGED KETLIB DATASET 
COMMON/KEYL I 3/KEYS (3258 1 

C0MM0N/LA0SAM/KARD(3),LIST(3),LA0SFC(3,6 l.DDNAMEdO) 
C0MH0N/TYPrS/03JNMF(24) 

C CQMMON/WTRKIN/SPACE t6*K0HPS*KG HPS+8*KOHPS | 

CQMMON/wr'RKIN/SPACFISaO) 

C OBJNME=THF ACTUAL OBJECT NiMElS) OF THE TYPE-OBJ ECT ( S ) TO B= PROCESS'^D 
C HR03J=NIJMR=o OF TYPE-OBJECT(S) TQ BE STATISTICALLY PROCESSED 
C 0BJFC=THE actual OBJSCT-NUNBERIS ) OF THF TYPE-OBJECTC S) IN THE TYPE-LIBRARY 
data CBJNME/'CLSARCTa* . ’CLEARC72 ' REGROWll • , 'REGR0WL6 't 'PINE » 

2/ 

DATA NPOBJ,OF1J=C/3,1 ,2.3f A,5/ 

DATA FPS, ITET8/1. 0-08*20/ 

DATA OnMAMF/i KEYLIB • , • OAT ASFT I */ 

DATA DELETE/'PINE '/ 

FND 


COMMON/WORKIN/ SPACE (x) 

DATA 0BJNME/'CLEARC73’ 

,’PINB V 

DATA NROBJ,OBJEC/5,l,2;3,4j5/ 


X must be at least (6*K0MPS*K0MPS+ 

+8*K0MPS) 

OBJNME(l), Names of the type- 

objects' to be processed, 

NROBJ Number of type-objects to 
be processed, 

OBJEC Actual object-Tiumber (s) 

(increasing values) of the 
type-objects in the type- 
library, 

OBJNME(l), must correspond to 
OBJEC(l), and be in 
numerical order. 
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DATA EPS,ITER8/1.D'08,20/ EPS Criterion 'of interruption. 

ITERS Number of iterations in 

eigenproblem calculations* 

DATA DELETE/^PINE V If ’PINE ’ exists as a type- 

class in KEYLXE, it will be dele- 
ted before any possible further 
processing is started. 

Subroutines in this task: 

INTRAK ( raTBYT,KEYXAP (LNDET , BYTINT , ESC ALB , FLOBYT , INTBYT , SCALER) ) 





6<36-lNTRA<PY PAGE 


12 


TYPF OBJECT: PJNE 

INPUT TYPBLIB; N584£i22 

LINEl 

LINES 

LOCUS 1 

TYPE-OBJFCT 4PEiUS): 

1291 

4 

1483 

INTRACLASS MFAh» VcCTor; 

4 l,99687f5F + 01 

5 K234375F+01 

7 1.584375F+01 

iT’'te 

- 

4 

1453 


(C3MgR:,0r AGONAL parti: 
*2 


rwrcACL<iS« OISPEPSION MdTRfX 
5 

4 4 , 677734^-01 

5 2,607422F-01 9.130Q^9F-01 

7 a.A'SinaF-Ol 1.162109?-01 Iv40336E+OO 


LOCI OF MSSBANDS 

0034507000 
4 0004507000 

4 


rHF FrGeNPRPRLrM f^LCULAnriMS 01 SC aWTI ^^U? 0 fiFTER 20 fTERATrOMSaR ^HEt4 |F( < 

genvcCtc^^^ np THE intpa:lass dispersion MAT*^IX: 

4 2, 749230F-0 1-2 • 955nZl F-0 i- 9» i 491 OOF-0 I 

5 2, 79L733F-0 1-6.S60414<=-01 3.701254P-01 

7 9,200432E-01 3.b71744S-01 1*610001F-01 


l*00-3tt IN the search FOR A ROOT OF THE eiGENFUNCT ION F 


^l*’'uooooooDvoo^°'''‘''*^'"’'* aowe» diagonal PART): 

5 -1.22120^0-10 1.0000000+00 
7 -2.79100GD-11-6.P75125D-U 1.0000000+00 


a\ 

o^ 


EIGENVALU!=S CF TH= INT^ACLAFS rilSPcsSTON 
1. ssas-fO'+oo 9.532i7i=-oi o.ionso^-oi 


MATRIX: 


-7IA99L1 =-01 f'F’'F»MlNANT OF HF INTRACLASS DISPERSION MATRIX: 


MAXIMUM LIKfLtH'-'On TLASSIFISP niTIAL TRAINING: 'VALUATION OF TVPE-CLASS STATISTICS QJANTITIES FOi? THE N684E122 KEY LIBRARY 
TH= NUMBER OF R EPO ES ENTAT T VES EXTRACTED FROM THE SCENE OBJECT TYPIFYING THIS CLASS I S 32 
THE 1M5S4EI22 KEyLIB now contains 5 TYPE-CLASSES WITH A TOTAL Of 140 REPRES ENT ATI VES 


NOW CONTAIflJ THE FOLLOWING TYPE-CLASSES 

2: CLEARC72 
3: RESPOHU 
4: REGR0H16 
5: PINE 

KROPPEFJALL, OALSLAND, clear-cut inventory (OPHAUGfWASTENSON) 
DATE 76-09-14 



67 


INTERKEY 


Evaluates final interobject statistics quantities to be used for 
ma^cinium likelihood classification. 


BLOCK DATA 

C FOa USE WITH F0A355 TELEGNOSTICS SUBROUTINE LIBRARY 
C 

REAL*8 DSNAHEfODNAMEtEPS 

COMMCM/TSLIB/DSNAME,UINEl,LINeS,L0CUSl,L0CItNSSB(10) 
COHMON/LA!5S«M/KARO(3l f L IST( 3 ) ,LAQSFC( 3f 6) ,DONAHE 1 10) 

CQMMON/eSCALA/FPSf ITERS 
f 

C C0MH0N/K5YLIB/KEYS( LENGTH) 

C LEWCTH MUST NOT BE LESS THAN THE BLOCKSIZE OF THE CATALOGED KEYLI3 DATASET 
COK«0N/KEYLlB/KEYSt3258) 

C CnMMnN/MORKIN/SPACE(6*KOMPS*KOHPS+8*K3HPS) 

COMHTNZ WORK IN/SPACE (680) 

DATA DDNAME/ 'KcYLIB »/ 

DATA EFS,ITER8/1.0-08flS/ 

FNO 


DATA EPS,ITER8/1.D-08,15/ EPS Criterion of interruption 

ITERS Number of iterations in 

eigenproblem calculations. 

Subroutines in this task: 

INTERK ( INTBYT , AVECT , DTOTAL , ESCAL8 , MATMUL , OMITRA , SECMOM) 


ORIGINAL PAGE IB 
OE POOR QUALITY. 
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TASK DOCUMENTATION LOG FOR FOA355 MULTILAyER DIGITAL IMAGE DATA HANDLING AND PROCESSING ARCHIVE (76-09-lA) 

MAXIMUM LIKELIHOOD CLASSIFIER FINAL TRAINING: EVALUATION OF InTERCLASS STATISTICS aUANTlTlES FOR THE N584E122 KEY LIBRARY 

INTERCLASS MEAN VECTOR: 

A 2.255713E+01 
5 I.50SO00F+01 

7 2.029997E+01 

INTERCLASS DISPERSION MATRIX (LOWER DIAGONAL PART}: 

A 2.975560F*00 
5 2.529A5T'^+00 A.40A651E+00 

7 A.2A7396E«-0O l.iA2282F+00 1.S73882E+-01 

THE EIGENPROBLEM CALCULATIONS DISCONTINUED AFTER 15 ITERATIONS OR WHEN iFl < l.DD-08 IN THE SEARCH FOR A ROOT OF. THE EIGENFUNCTION F 

EIGENVECTOR MULTIPLICATION CHEC< = UNIT MATRIX (LOWER DIAGONAL PART): 

A l.OOOOOOD+00 
5 -3.5704170-13 UOOOOOOD+00 
7 -4.0280740-13-3.0949960-12 1.300000000 

EIGENVALUES OF THE TRANSFORMFD OM IT-CLEARC73 DISPERSION MATRIX; 

2.212044F+01 4 .771733EFOO 1 .1 58403FFOO 

EIGENVECTOR MUL TI PLI CAT! 3N CH>^C< « UNIT MATRIX (LOWER DIAGONAL PART): 

4 l.OOOOOOD+DO 

5 -1.07706BD-11 1 .00000004-00 

7 -1.3877790-16 1.2528230-11 1.0000000+00 

EIGENVALUES OF TH^ TRANSFORMED OM IT-CLEARC72 DISPERSION MATRIX: 

2.046512E+01 6.578749E+00 5.212763F-01 

EIGENVECTOR MULTIPLICATION CH^CX = UNIT MATRIX (LOWER DIAGONAL PART): 

A 1. 0000000+00 

5 3.759877D-13 l.OOOOOOD+00 

7 -A. 13892AD-14-1.09.528A0-11 l.OOOOOOD+00 

EIGENVALUES OF THF TRANSFORMED OMlT-REGROWll DISPERSION MATRIX: 

2.0930A7E+01 1.00375AE+01 8.9a5035E~01 

'j ‘ ^ 

EIGENVECTOR MULT IPL TCA TIDN CH'^CX = UNIT MATRIX (LOWER DIAGONAL PART): 

4 1 .0000000+00* > 

5 1. 196265D-M‘r.OO0O00D+0O 

7 -1. 124101D-1 5-1. 6031620-13 1. 0000000+00 

EIGENVALUES OF-*,THE TRANSFORMED DM IT-REGRDW16 DISPERSION MATRIX: 

1.0681S1E+01‘2.362064E+00 5-911416E-01 

EIGENVECTOR MULTr'iPLICATION CHEC< = UNIT MATRIX (LOWER DIAGONAL PART): 

4 1.0000000.+.00/ 

5 -6.58248TD-12£.1.000000D+00 

7 -4.830336D-11-1.011840D-10 l.OOOOOOD+OD 


EIGENVALUES OF THE TRANSFORMED OMIT-PINE -DISPERSION MATRIX: 

4.4A0283E+01 9.164378E+00 1.792923E+00 

the N584E122 KEYLI6 CONTAINS 5 TYPECLASSES WITH APPROPRIATE STATISTICS QUANTITIES EVALUATED AND READY FOB USE IN CLASSIFICATIONS 
KROPPEFJALL, OALSLANO, “cLEAR-C’uT inventory (ORHAUGtWASTENSON) 
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TYPESORT 


Maximum likelihood classification of a selected window using t;he 
statistically preprocessed data of the appropriate KEYLTB* 


BLOCK DATA 

C POR USF WITH FOA355 TELEGNOSTI CS’ SOBROOTIH E LIBRART 

C 

*"*?EAL’^8 DSr^AMFfDDNAME 

COMMON/TGLI8/DSNAHE,LINEltLIMEStLOCUSULOCI,HSS3t 10 \ 
COHN^nN/LAOSAM/KAROO) tLISTO) ,LAQSFCt 3,6 l,DDWAMEI 10> 
COM-^ON/MAKOev/BDUNO 

C 

C COMMON/KEYLIB/KEYSaEMCTH) 

C UEMGTH MUST NOT BE LESS TH&M THE BLOC<SIZE OF THE CATALOGED KEYL13 DATASET 
COMMON/ KEYL I BTKEYS (3258) 

C rnMM0N/W0RKTN/SPACE(6*K‘)MPS*K0MOS<-8*K3MPS) 

COMMON/ WOP< IN/S PACE(680) 

C 

DATA DSM&ME/ ‘SORTED V 
DATA BOUND/2 *0/ 

DATA LINEl, LINES, LHCUSl, LOCI ,MSSB/ 120 1 ,l 28 ,1 1 123, 4,5, 7/ 

DATA DDNAHF/‘KFYLI6 S ‘ DATASET I •,» SORTED •/ 

END 


COMMON/WORKIN/ SPACE (x) x must be at least (6*KCfMPS*K0MPS+ 

+8*KOMPS) 

DATA BOUND/2.. 0/ Using 2,0 standard deviations per 

channel in the classification as 
acceptance/rejection criterion, 

DATA LINEl,..,. See'the description of the source 

deck. 

Subroutines in this task: 

MSSORT (INTBYT, CULLER(TGHEAD,INTBYT,FLOBYT, BYTINT)) . 
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^8fi-TYPESnPT PAGB^ 

TASK OOCUMPMTATI0fi i*QG PHR ^0*355 MULTILAYER DIGITAL IMAGE DATA HANDLING AMO PROCESSING ARCHIVE {76-09-lA) 
LANDSAT: 4 


scene/frame id 

^''EXPQSUPE 

SCFNH CENTER 

ORIGINAL COftlQLlSRECTIFIED 


sodd-hhmmspm 

*o‘ate; jt 

L4T LONG 

CCTCS): 

LINE LENGTH 


1330-0Q523Q0 

18JUN73 C 

M58-41/E012- 

23 1 

2 3 A 

3358 


INPUT DATfiSET 

SUBSCcME: 

LlNEl 

LINES . 

LOCUSl 

< LOCI OF 

MSSBANDS 

K^OPPTFJ 


1073 

256 

13A3 

256 A 

5 6 7 

CLASSTFtRO SUOSCEME; SORTFO 1201 

128 

lAAA 

128 

ALL 3 


MAXIMUM LIKEirHOOD CLASSI F IC ATION USING THE N584E122 KEY LIBRARY 

sorting into 5 CLASSES BUT REJECTING EVERY PIXEL OUTSIOE 2*0 STANDARD OEVlATIONtS) PER HSSBAMO FROM ITS MAXIMUM LIKELIHODO CLASS 
0! REJECTS : 

15 CLEARC73 
2: CLEARC72 
3: REG^OWli 
A: REGROWlf 
5; PINE 

KROPPEEJilL, DALSLANl), CLEAt-CUT INVENTORY { DRHAUGi WA STENSON) 


DATE 76-09-14 
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GROUPINP JF TASKS BY I/O PMITS 



CCTMERGE 

REPLIC8 

SIXPACK 



MERGTAP 

CCTAP 

CNESTAP 



NUANCES 
HISTOGR 
FOA 355 
FOA 356 
COVCORR 
HEXPIXEL 



CALCOMP (microfilm plotter) 


PROJEC2D 

CALCOMP 


ORIGINAL PAGE 
OF POOR QUALITY 
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Hertz (ink plotter) 



Ronntek (TV-monitor) 


INKJETS 

RAMTEK 



COMPOUND 

MOVE 

TYPEHEAD 

TYPELIB 

KEYHEAD 

INTRAKEY 

INTERKEY 

TYPESORT 
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TYPICAL APPLICATION TASKS 
IMAGE GENERATION 



OP SOOR 
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TRAINING AND CLASSIFICATION 

(1t^ 



FOA3S6 
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0BUjr-5JAL page is 
OF Poor quality 




















